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I.  INTRODUCTION 


For  about  the  past  six  years,  high-current  betatrons  of  various 
designs  have  been  the  subject  of  experimental  and  theoretical  Investiga¬ 
tions.  This  Interest  stems  from  the  potential  for  compact,  high-power  beam 
generation  which  these  machines  offer. 

To  operate  a  high-current  betatron  successfully,  one  must  be  able  to 
inject  the  beam  onto  a  closed,  recirculating  orbit,  and  maintain  beam 
stability  for  the  duration  of  the  acceleration  cycle.  The  latter  may  vary 
from  tens  of  microseconds  to  several  milliseconds.  Depending  on  the  appli¬ 
cation,  it  may  then  be  necessary  to  extract  the  beam  from  the  device. 

During  the  period  of  our  contract  with  the  Office  of  Naval  Research 

(01  November  1983  -  31  March  1987),  we  have  concentrated  on  the  question  of 

beam  behavior  after  injection,  because  this  behavior  determines  such  key 

parameters  as  the  maximum  operating  current,  and  the  length  of  the  acceler¬ 
ation  cycle.  During  the  period  01  November  1983  -  31  December  1984,  we 
developed  the  analytic  and  computational  tools  to  study  the  equilibrium  and 
stability  of  the  circulating  electron  ring.  This  work  was  described  in 
detail  in  Ref.  1,  and  will  be  referred  to  as  needed.  In  the  period 
01  January  1985  -  31  March  1987,  which  Is  covered  in  this  report,  our 
emphasis  was  on  applying  these  methods  to  experimental  machines  currently 
in  operation. 

The  experiments  for  which  we  have  performed  specific  calculations  are 
the  Modified  Betatron  Accelerator  at  the  Naval  Research  Laboratory  [2],  the 
Stellatron  at  UC  Irvine  [3],  and  the  Solenoldal  Lens  Betatron  at  the 
University  of  New  Mexico  [4],  This  report  is  organized  Into  three  main 
sections  describing  the  results  for  each  of  these  devices.  Publications, 
Conference  Proceedings  and  Presentations,  and  Technical  Reports  written 
during  this  contract  are  listed  in  Appendix  D. 


II.  MODIFIED  BETATRON  ACCELERATOR 


1.  NEGATIVE -MASS  INSTABILITY  ON  A  COLD  BEAM 

The  modified  betatron  adds  a  toroidal  magnetic  field  to  the  conven¬ 
tional  betatron  configuration.  In  this  way,  the  space-charge  limit  on  the 
confinable  current  is  considerably  increased.  The  most  important  collec¬ 
tive  instability  for  this  device  is  the  negati ve-mass  instability,  which 
causes  toroidal  clumping  and  transverse  kinking  of  the  beam.  In  Ref.  1,  we 
derived  a  dispersion  relation  for  this  Instability  which  agrees  accurately 
with  cold-beam  numerical  simulations.  In  Fig.  1,  this  dispersion  relation 
is  applied  to  the  parameters  of  the  modified  betatron  at  NRL.  The  present 
status  of  this  machine  is  that  a  1  kA,  0.8  MeV  beam  has  been  successfully 
trapped  in  the  combined  vertical  and  toroidal  magnetic  fields  for  periods 
up  to  10  usee.  During  this  time,  the  beam  executes  on  the  order  of  500 
major  revolutions.  The  observed  stability  is  in  agreement  with  the  predic¬ 
tion  of  Fig.  1.  Due  to  space-charge  depression  of  the  Injected  beam 
energy,  the  actual  value  of  the  relativistic  factor  y  for  the  circulating 
ring  Is  about  1.6  [5].  This  puts  the  beam  below  the  "negative-mass 
transition  energy",  which  is  given  by 


where  v  is  Budker's  parameter,  R  Is  the  major  radius,  and  a  is  the  torus 
minor  radius.  This  gives  y^  •  2.1  for  a  1  kA  beam,  in  agreement  with 
Fig.  1.  In  the  absence  of  stabilizing  measures,  the  beam  will  become 
unstable  above  this  energy.  The  peak  growth  for  the  i  ■  1  instability 
occurs  around  y  ■  4,  and  the  mode  would  become  nonlinear  in  a  time  on  the 
order  of  1  usee.  We  have  not  simulated  the  nonlinear  development  of  the 
negative-mass  Instability  for  a  1  kA  beam.  For  10  kA  beams,  our  simula¬ 
tions  show  strong  kinking  of  the  beam  at  low  energies  leading  to  considera¬ 
ble  current  loss  (see  Fig.  4a  below,  and  Ref.  6). 


2 


BETATRON  ' 

v 


.  Growth  rates  of  the  negative-mass  instability  for  parameters 
typical  of  the  NRL  betatron. 


2.  EFFECT  OF  SPREAD  IN  CIRCULATION  FREQUENCIES 


One  proposal  to  reduce  or  eliminate  the  negative-mass  instability 
involves  a  spread  in  particle  circulation  frequencies  [7].  Some  spread  in 
circulation  frequencies  is  always  present  due  to  finite  beam  radius,  but 
since  the  beam  rotates  poloi dally  [8],  the  frequency  spread  is  only  of 

p 

order  (rb/R)  ,  where  rb  is  the  beam  radius.  For  a  given  beam  radius,  one 
obtains  maximum  frequency  spread  when  the  beam  particles  do  not  oscillate 
transversely,  as  discussed  in  Ref.  1.  One  can  achieve  this  by  introducing 
an  energy  spread  on  the  beam,  or  by  making  the  poloi dal  rotation  frequency 
small  enough.  The  latter  occurs  naturally  near  the  so-called  diamagnetic- 
paramagnetic  transition  energy,  where  the  poloidal  rotation  changes  direc¬ 
tion  [8]. 


For  the  low-current  conventional  betatron,  it  is  possible  to  derive 
analytically  the  effect  of  a  spread  in  circulation  frequencies.  For  the 
case  of  a  beam  with  a  uniform  current  density  profile,  one  finds  [see 
Appendix  C,  Eqs.  ( 34 ) - ( 38) ]  that  the  growth  rate  goes  to  zero  for  a  beam 
radius  such  that 


c  >  r. 


(2) 


where  rQ  is  the  growth  rate  for  a  cold  beam  and  t  is  the  toroidal  mode 
number.  If  the  beam  has  a  finite  poloidal  rotation  frequency  “p0^ »  Then 
one  can  argue  heuristically  that  Eq.  (2)  is  still  the  condition  for  sta¬ 
bility  provided,  that,  in  addition 


This  extra  condition  ensures  that  a  disk  of  particles  initially  at  one 
toroidal  location  will  become  smeared  out  over  several  wavelengths  during 
one  poloidal  rotation  period. 
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We  have  performed  numerical  simulations  to  test  the  applicability  of 
;  Eqs.  (2)  and  (3)  to  the  modified  betatron,  for  which  we  have  not  attempted 

a  rigorous  derivation  of  these  criteria.  We  looked  at  cases  in  which  the 
beam  current  and  energy  are  kept  fixed,  but  where  the  poloidal  rotation 
frequency  of  the  beam  is  varied.  This  is  accomplished  by  changing  the  beam 
radius.  The  two  sides  of  Eq.  (3)  are  plotted  versus  beam  radius  in  Fig.  2 
for  a  10  kA,  y  =  12  beam.  For  a  beam  radius  of  2  cm,  Eq.  (3)  is  strongly 
violated  for  low  t  modes  (we  have  chosen  t  =  4  for  this  example).  Thermal 
effects  are  therefore  expected  to  be  negligible.  This  is  indeed  what  we 
find  in  the  simulation.  As  shown  in  Fig.  3a,  the  instability  grows  at  the 
predicted  cold  growth  rate.  A  particle  plot  made  near  the  end  of  the 
linear  growth  is  shown  in  Fig.  4a.  Nine  disks  of  particles  are  used  to 
resolve  one  wavelength  of  the  instability.  Initially,  these  disks  have 
zero  thickness,  but  a  spread  in  circulation  frequencies  causes  them  to 
smear  out  toroidal ly.  Consistent  with  the  observed  cold  beam  growth  rate, 
the  disks  in  Fig.  4a  show  little  spreading.  Soon  after  the  stage  in 
Fig.  4a,  about  3/4  of  the  current  is  lost  to  the  walls.  Parenthetically, 
we  remark  that  this  simulation  illustrates  that  the  toroidal  magnetic  field 
is  a  mixed  blessing  for  beam  stability.  While  it  cuts  the  growth  rate  of 
the  instability  significantly  from  what  one  would  obtain  without  it,  it 
also  inhibits  the  radial  motion  which  tends  to  saturate  the  instability  by 
(  producing  a  spread  in  circulation  frequencies.  Note  that  the  distortion  of 

the  beam  in  Fig.  4a  is  primarily  in  the  z  direction,  and  so  does  not 
contribute  to  stability. 

In  contrast  to  the  case  just  discussed,  Eq.  (3)  is  well  satisfied  when 
the  beam  radius  is  increased  to  3.5  cm,  so  that  the  poloidal  rotation  may 
be  effectively  ignored  (see  Fig.  2).  Substituting  the  cold  beam  growth 
rate  into  Eq.  (2),  we  find  that  marginal  stability  is  predicted.  In  a 
simulation  for  this  beam  radius,  however,  we  see  a  finite  growth  rate  of 
about  half  the  cold  growth  rate,  as  shown  in  Fig.  3b.  Thus,  Eq.  (2)  some¬ 
what  underestimates  the  beam  radius  required  for  stability  in  the  modified 
betatron.  The  instability  causes  the  beam  to  distort  is  a  manner  similar 
to  the  small -radius  case,  as  shown  in  Fig.  4b.  However,  the  instability 
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Figure  2.  Plot  of  the  two  sides  of  Eq.  (3)  for  a  10  kA,  y  =  12  beam 
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Figure  3.  Energy  histories  of  the  £.  =  4  negative-mass  instability  on  a 
10  kA,  y  =  12  beam  with  initial  radius  r^  =  2  cm  in  (a)  and 
rj,  *  3.5  cm  in  (b). 
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Figure  4,  Particle  plots  taken  during  the  nonlinear  stages  of  the  Simula 
tlons  shown  In  Fig.  3.  For  (a),  the  Initial  beam  radius  Is 
rK  *  2  cm  and  In  (b),  the  Initial  beam  radius  Is  rK  ■  3.5  cm. 


saturates  at  a  much  lower  level  for  the  larger  beam  radius.  This  leads  us 
to  believe  that  the  beam  is  close  to  marginal  stability  initially,  and  that 
the  small  energy  spread  created  by  the  instability  is  sufficient  to 
stabilize  it  without  loss  of  current. 

We  have  run  a  similar  series  of  cases  for  the  I  =  20  mode.  These  runs 

illustrate  again  the  importance  of  satisfying  Eq.  (3)  in  order  to  get  the 

full  benefit  of  finite  beam  radius  stabilization.  For  a  1.7  cm  beam 

radius,  we  see  from  Fig.  2  that  Eq.  (3)  is  not  satisfied  for  the  t  =  20 

mode.  A  particle  simulation  shown  in  Fig.  5a  confirms  that  the  instability 

-3  -1 

grows  at  the  cold  beam  growth  rate,  Tg  *  2.9  x  10  cm  .  Increasing  the 
beam  radius  to  2.0  cm  puts  us  in  a  regime  where  Eq.  (3)  is  only  weakly 
satisfied.  We  observe  growth  of  the  instability  at  roughly  half  of  the 
cold  value,  as  shown  in  Fig.  5b.  The  wavy  energy  curve  seems  to  be  charac¬ 
teristic  of  cases  where  Eq.  (3)  is  weakly  satisfied.  Note  that  Eq.  (2) 
predicts  stability  for  the  2  cm  beam  (irb/R  =  4  x  10  cm  ).  The  non¬ 
linear  development  is  very  similar  to  that  of  the  1.7  cm  beam. 

Increasing  the  beam  radius  further  to  3  cm  results  in  Eq.  (3)  being 

well  satisfied,  and  Eq.  (2)  becomes  satisfied  by  a  factor  of  two  ( rnc  = 

-3  -1  2-3-1  u 

2.9  *  10  cm  ,  trfe/R  =  6  x  10  cm  ).  In  a  simulation  for  this  radius, 

some  erratic  growth  is  seen  (Fig.  5b),  but  not  enough  to  measure  a  reliable 

growth  rate.  The  peak  perturbed  field  energy  in  Fig.  5b  is  only  0.2%  of 

the  equilibrium  field  energy.  The  beam  radius  increases  slightly,  but  no 

current  loss  is  observed. 

Finally,  we  present  some  results  for  parameters  closer  the  the  NRL 
experiment  than  the  above  cases.  We  performed  a  simulation  of  the  l  -  20 
mode  near  the  peak  of  the  growth  curve  In  Fig.  1  for  a  3  kA,  y  *  7  beam. 

The  two  sides  of  Eq.  (3)  are  plotted  in  Fig.  6a.  We  choose  a  radius  of 
2  cm, which  is  on  the  order  of  the  experimentally  measured  value.  For  this 
radius,  Eq.  (3)  is  marginally  satisfied,  and  Eq.  (2)  is  satisfied  by  about 
a  factor  of  three  (rQc  »  1.3  x  10"3  cm"*,  Ar^/R^  *  4.0  x  10"3  cm"1).  The 
growth  rate  observed  in  the  simulation  is  erratic,  and  on  the  order  of 
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Figure  5.  Energy  histories  and  particle  plots  for  simulations  of  the 

i  *  20  negative  mass  Instability  on  a  10  kA,  y  »  12  beam.  The 
Initial  beam  radius  Is  1.7  cm  In  (a),  2  cm  In  (b)  and  3  cm  In 
(c). 
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Figure  6.  Plot  of  the  two  sides  of  Eq.  (3)  for  parameters  typical  of  the 
NRL  betatron.  In  (a)  the  beam  parameters  are  3  kA,  y  =  7,  and 
in  (b)  the  parameters  are  1  kA  and  y  •  4.  In  each  case,  a 
2 kG  toroidal  field  is  assumed. 
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3.6  x  10  cm  ,  much  less  than  the  cold  growth  rate.  If  Eq.  (3)  were 
better  satisfied,  e.g.,  rfe  _>  3  cm,  the  t  =  20  mode  would  likely  be 
completely  stable. 

In  Fig.  6b,  we  show  the  two  sides  of  Eq.  (3)  for  a  1  kA  beam  in  a 
2.5  kG  toroidal  field  at  an  energy  of  y  a  4,  which  corresponds  roughly  to 
the  peak  of  the  1  kA  growth  rate  curve  In  Fig.  1.  We  see  that  for  low 
t -numbers,  Eq.  (3)  is  not  satisfied  for  rb  £  2  cm.  We  therefore  would 
expect  to  see  some  growth  of  the  negative-mass  instability  at  these  mode- 
numbers  in  the  NRL  experiment.  The  nonlinear  development  of  the  insta¬ 
bility  is  not  known  for  a  1  kA  beam  current.  Given  the  large  minor  radius 
of  the  drift-tube  (15  cm).  It  Is  possible  that  the  instability  may  saturate 
without  loss  of  current. 

3.  BEAM  STABILITY  AT  HIGH  ENERGIES 

A  novel  prediction  of  the  analytic  dispersion  relation  derived  in 
Ref.  1  is  that  the  character  of  the  negative  mass  instability  for  high  i 
numbers  changes  above  a  certain  energy.  "High  l"  in  this  context  typically 
means  l  >  12  [1].  A  representative  growth  curve  is  shown  in  Fig.  7  (i  = 

20,  I  =  10  kA,  Be  =  1  kG).  Above  y  *  20,  the  usual  nonresonant  negative 
mass  instability  gives  way  to  a  resonant  Interaction  between  the  longi¬ 
tudinal  and  transverse  modes  of  the  beam.  We  performed  a  number  of  simula¬ 
tions  in  this  regime  to  test  the  analytic  predictions.  The  results,  shown 
in  Fig.  7,  are  in  reasonably  good  agreement  with  theory.  A  common  feature 
of  the  simulations  for  y  >  20  is  that  the  beam  radius  increases  in  the  non¬ 
linear  state,  but  no  current  loss  occurs.  The  nonlinear  state  of  the  beam 
in  the  y  *  50  simulation  Is  shown  In  Fig.  8.  A  steady  state  transverse 
profile  is  not  reached  In  the  time  of  the  simulation.  However,  the  energy 
In  the  perturbed  fields  decreases  by  about  a  factor  of  30  from  its  peak 
value  during  the  nonlinear  saturation,  as  the  bunches  of  charge  become 
smeared  out.  The  relative  ease  with  which  the  modes  stabilize  nonlinearly 
can  be  attributed  In  part  to  their  relatively  low  growth  rates,  and  to  the 
fact  that  little  transverse  kinking  of  the  beam  is  produced  by  the 
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Figure  7.  Instability  growth  rates  of  the  l  *  20  mode  for  a  10  kA  beam  In 
a  1  kG  toroidal  field. 
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Instability,  unlike  the  negative-mass  instability.  This  is  discussed 
further  below. 

For  the  simulations  at  y  a  35  and  50,  the  two  sides  of  Eq.  (3)  are 
comparable  in  magnitude,  so  that  it  is  difficult  to  predict  how  much 
damping  from  circulation  frequency  spread  is  present  in  the  simulations. 

If  Eq.  (3)  were  well  satisfied,  then  Eq.  (2)  would  predict  stability  for 
each  of  the  simulations  for  y  >  20.  Based  on  the  results  in  Figs.  4b  and 
5b,  it  is  perhaps  surprising  that  there  is  not  a  larger  discrepancy  between 
the  cold  beam  growth  rates  and  the  simulation  results. 

For  the  simulation  at  y  =  20,  Eq.  (3)  is  reasonably  well  satisfied,  so 

that  we  would  certainly  have  expected  the  thermal  effects  to  have  a  more 

significant  effect  than  is  evident  from  Fig.  7.  Since  the  shape  of  the 

analytic  curve  in  this  region  is  sensitive  to  small  changes  in  the  beam 

energy,  we  also  performed  simulations  at  y  =  19.7,  and  y  =  20.5,  but 
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obtained  the  same  growth  rate  as  for  y  *  20,  namely,  rQ  =  1.1  x  10  cm  . 

Next,  we  sought  to  decrease  whatever  thermal  effects  might  be  present  by 

reducing  the  beam  radius  from  2  cm  to  1  cm.  This  made  the  beam  rotate 

faster  In  the  pololdal  direction.  The  growth  curve  for  this  beam  at  y  =  20 

is  shown  in  Fig.  9.  Note  that  there  are  two  distinct  slopes,  giving  growth 
-4  -1  -3  -1 

rates  of  6.3  *  10  cm  and  1.56  x  10  cm  .  It  Is  tempting  to  identify 
the  lower  growth  rate  with  the  dip  in  growth  rate  at  y  =  20  in  Fig.  7.  It 
may  be  that  the  instability  grows  at  the  analytically  predicted  rate  at  low 
amplitudes,  and  that  a  type  of  mode-locking  takes  over  as  the  amplitude 
gets  large.  In  this  light.  It  may  be  significant  that  in  each  of  the  2  cm 
beam  runs  around  y  *  20  there  Is  a  relatively  long  period  (-100  ns)  at  the 
start  of  the  simulations  where  very  little  growth  is  seen.  At  Jt  =  20,  the 
resonant  denominator  In  the  dispersion  relation  due  to  the  TE^  mode  is 
starting  to  have  a  noticeable  effect  on  the  beam  behavior,  even  though  we 
are  not  at  exact  resonance.  This  Is  evident  from  the  TE^-like  structure 
of  the  contour  plot  of  the  perturbed  magnetic  field  B.  shown  in  Fig.  10b, 
which  Is  taken  from  the  nonlinear  stage  of  the  y  *  20,  r^  =  1  cm  simula¬ 
tion.  It  is  Indeed  the  proximity  of  the  TE^  mode  resonance  which  is 
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responsible  for  the  change  in  the  character  of  the  negative-mass  insta¬ 
bility  referred  to  above. 

For  the  1  cm  beam  case,  there  is  strong  longitudinal  bunching  of  the 
beam  in  the  nonlinear  state,  as  shown  in  Fig.  10a.  This  contrasts  with  the 
2  cm  beam  runs,  where  the  perturbations  saturate  at  a  much  lower  level.  It 
also  contrasts  with  the  nonlinear  development  at  y  =  12,  shown  in  Fig.  4a, 
where  there  is  considerable  transverse  kinking  of  the  beam.  Such  pro¬ 
nounced  transverse  kinking  is  characteristic  of  the  negative-mass  insta¬ 
bility.  The  different  nonlinear  behavior  is  evidence  that  the  instability 
for  y  >  20  is  not  a  negative-mass  instability.  It  is  possible  to  have 
longitudinal  bunching  without  much  transverse  kinking  since  the  equilibrium 
position  of  the  beam  depends  on  y  +  where  4  is  the  electrostatic  poten¬ 
tial  of  the  beam;  y  +  can  remain  fixed  as  the  beam  bunches. 

We  also  ran  a  3  cm  radius  beam  at  y  =  20.  Equations  (2)  and  (3) 
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predict  stability,  but  we  observe  a  growth  rate  of  8.2  x  10  cm  ,  only 
slightly  smaller  than  the  2  cm  rate.  The  mode  saturates  at  approximately 
the  same  level  as  the  2  cm  runs.  Finally,  we  ran  a  2  cm  beam  case  where, 
instead  of  a  slow  E  x  B  mode  equilibrium,  we  ran  the  beam  in  the  Brillouin 
limit.  In  this  limit,  the  beam  is  rotating  at  the  Larmor  frequency,  0./2y, 
and  the  particles  are  given  whatever  temperature  is  required  to  give  force 
balance.  This  results  in  a  beam  with  little  spread  in  particle  frequen¬ 
cies.  We  observe  a  low  initial  growth  rate  of  4.2  x  10”4  cm"1,  (cf.  the 
1  cm  radius  case  described  above)  and  the  perturbations  saturate  at  an 
energy  level  about  an  order  of  magnitude  lower  than  for  the  slow  E  x  b 
equilibrium.  A  comparable  increase  in  beam  radius  is  observed,  however, 
whether  due  to  the  Instability,  or  to  slow  (perhaps  numerical  due  to  the 
fast  rotation)  changes  In  the  beam  equilibrium. 

4.  RESONANCES  IN  THE  MODIFIED  BETATRON  AT  LOW  ENERGIES 

In  addition  to  the  collective  instabilities  discussed  above,  the 
modified  betatron  is  prone  to  single-particle  resonances.  These  occur  when 
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one  of  the  transverse  oscillation  frequencies  is  an  integer  multiple  of  the 
circulation  frequency,  u>g: 

n±  -  t«g  (4) 


Errors  in  the  applied  magnetic  fields  then  result  in  transfer  of  energy 
from  the  longitudinal  direction  to  the  transverse  direction.  The 
transverse  frequencies  are  given  by  the  expressions 
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where  u>  is  the  focusing  force  due  to  betatron  focusing  and  self-fields. 

r  2  2  2  2 

For  the  beam  centroid  motion,  is  given  by  (1/2  -  n  r‘/a V,  where  n  is 

2  3  2  s  d  u  s 

the  self-field  index  (■>  2vR  /y  r£).  At  low  energies,  the  low  frequency  n_ 

is  less  than  u>Q,  so  that  resonance  is  impossible,  while  the  fast  frequency 

is  much  greater  than  u>n,  giving  a  weak,  high  ?rder  resonance.  For  the 

u  2  2 

motion  of  particles  within  the  beam,  on  the  other  hand,  ay  =  (1/2  -  ns)uQ, 

and  the  slow  frequency  can  be  a  few  times  larger  than  u>Q,  giving  rise  to  a 
low  order  resonance.  (The  fast  frequency  is  still  much  larger  than  o^.) 

The  ratio  |ft.|/oj0  is  plotted  as  a  function  of  energy  in  Fig.  11  for  param¬ 
eters  relevant  to  the  NRL  betatron  (1-3  kA,  BQ  =  2.5  kG,  rb  =  1  cm).  Near 
the  present  injection  parameters,  which  yield  a  1  kA  ring  with  y  =  1. 5-2.0, 
there  are  a  number  of  low  order  resonances.  Their  effect  on  the  beam  seems 
to  be  rather  weak,  however.  We  show  in  Fig.  12  the  effect  of  errors  of  1% 
and  4%  in  the  vertical  field  on  a  particle  near  the  I  -  2  resonance  in  a 
2  kA  beam.  The  particle  is  started  near  the  origin,  and  its  betatron 
motion  is  amplified  to  just  0.5  cm  for  the  4%  error.  It  is  not  possible 
for  the  motion  reach  an  amplitude  much  larger  than  the  beam  radius  since 
the  potential  well  due  to  the  self-fields  outside  the  beam  is  strongly 
anharmonic.  Also,  the  cases  shown  in  Fig.  12  are  not  self-consistent  in 
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Figure  12.  Amplification  of  betatron  oscillations  due  to  an  l  =  2  reso¬ 
nance  in  a  2  kA,  y  -  3  beam.  In  (a),  the  error  in  the  vertical 
field  is  1%,  and  in  (b)  the  error  is  4%. 


that  the  space-charge  of  the  beam  Is  held  fixed  while  the  test  particle  is 

being  advanced.  In  a  self-consistent  case,  the  whole  beam  would  expand, 

2 

thereby  changing  wj;  and  detuning  the  resonance.  We  therefore  believe  that 
these  resonances  may  contribute  to  some  heating  of  the  beam,  but  are 
unlikely  to  be  responsible  for  the  gradual  current  loss  observed  in  the 
experiment  [2], 
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III.  STELLATRON  ACCELERATOR 


1.  RESULTS  RELATED  TO  THE  UCI  STELLATRON 

The  stellatron  accelerator  adds  a  helical  quadrupole  to  the  modified 
betatron  design.  The  main  effect  of  this  is  to  increase  the  energy  accept¬ 
ance  of  the  device.  This  helps  to  suppress  the  negative-mass  instability 
by  inhibiting  transverse  motion  of  the  beam.  It  also  pushes  up  the  transi¬ 
tion  energy  at  which  the  effective  mass  of  the  particles  becomes  negative 
[9].  As  reported  in  Ref.  1,  our  simulation  results  support  these  con¬ 
clusions.  In  Fig.  13,  we  show  results  of  two  simulations  confirming  the 
increase  in  the  transition  energy  relative  to  the  modified  betatron.  The 
stellatron  growth  rates  in  this  figure  were  obtained  by  putting  the 
smoothed  approximation  for  the  quadrupole  focusing  obtained  by  Chernin  [9] 
into  the  modified  betatron  dispersion  relation  obtained  in  Ref.  1. 

The  UC  Irvine  stellatron  has  succeeded  in  accelerating  significant 
amounts  of  charge  from  rest  up  to  10  MeV  [10].  The  peak  current  is  over 
1  kA  up  to  about  4  MeV  and  gradually  drops  off,  leaving  a  few  hundred 
amperes  at  10  MeV. 

We  have  applied  the  dispersion  relation  derived  in  Ref.  1,  modified  as 
discussed  above,  to  the  parameters  of  this  experiment.  In  addition  to  the 
quadrupole  focusing,  however,  we  need  to  take  account  of  the  fact  that  the 
UCI  beam  is  charge-neutral.  This  is  because  it  Is  created  by  runaway 
acceleration  from  an  Initially  stationary  plasma.  The  runaway  fraction  is 
thought  to  be  on  the  order  of  100%  [10],  Thus  the  ions,  in  addition  to 
providing  space-charge  neutralization,  also  provide  transverse  focusing  of 
the  electron  beam,  in  the  manner  of  IFT  (ion-f.ocused  transport).  The 
transverse  focusing  frequency  becomes 

2  2,22 

«>r  ♦  *>r  +  y  nsa>0  (6) 
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Figure  13.  Growth  rates  of  the  l  *  20  negative-mass  instability  for 

typical  stellatron  parameters.  The  solid  line  is  for  a  stella- 
tron  field  such  that  the  helical  field  parameter  e  =  0.7,  while 
the  dashed  line  is  for  e  *  0. 
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When  we  include  the  additional  focusing  term  in  the  dispersion  relation,  we 
find  that  its  effect  on  the  beam  stability  is  greater  than  the  effect  of 
the  helical  quadrupole,  as  shown  in  Fig.  14.  In  Fig.  14a,  we  see  that 
introducing  the  ions  makes  a  significant  difference  in  the  peak  growth 
rate,  and  in  the  transition  energy.  In  Fig.  14b,  on  the  other  hand,  we  see 
•  that  applying  the  helical  quadrupole  has  little  effect.  (The  helical  field 

is  essential  to  generate  the  beam  Initially,  however.)  This  result  was 
suggested  by  Ishizuka  [10],  and  is  in  agreement  with  the  experimental 
findings.  Disruption  of  the  beam  accompanied  by  strong  rf  oscillations  at 
I  harmonics  of  the  circulation  frequency  is  seen  early  In  the  acceleration 

cycle.  This  is  thought  to  be  due  to  the  onset  of  the  negative  mass  insta¬ 
bility  as  the  beam  goes  through  the  transition  energy.  The  transition 
energy  is  observed  to  increase  with  the  beam  current,  which  is  consistent 
I  with  the  ion-focusing  picture.  For  a  1  kA  beam,  the  onset  of  instability 

is  observed  at  around  2  MeV,  in  agreement  with  the  prediction  from  Fig.  14. 

2.  ELECTROMAGNETIC  INSTABILITY  DUE  TO  THE  HELICAL  FIELD 

I 

A  virulent  electromagnetic  instability  observed  in  simulations  of  the 
stellatron  was  reported  In  Ref.  1.  We  subsequently  derived  an  analytic 
model  of  this  interaction  which  is  reproduced  in  Appendix  A.  The  mode  is  a 
»  three-wave  coupling  resembling  the  free-electron  laser  Interaction,  except 

that  the  beam  mode  involved  is  transverse  rather  than  longitudinal.  Using 
FEL  terminology,  the  quadrupole  field  field  acts  as  a  wiggler  and  upshifts 
the  slow  transverse  mode  on  the  beam  so  that  it  can  intersect  the  lowest 
#  electromagnetic  mode,  the  TE^  mode.  This  is  illustrated  qualitatively  in 

Fig.  15. 

The  analytic  model  given  In  Appendix  A  gives  reasonably  good  agreement 
I  with  the  simulation  results.  For  high  currents  and  strong  focusing,  the 

growth  rates  can  be  very  large,  A  simulation  of  a  10  kA,  3.3  MeV  beam 
showed  disruption  of  the  beam  In  just  one  turn  around  the  1  m  radius 
device.  We  have  seen  linear  growth  of  the  instability  In  simulations  with 
I  currents  as  low  as  300  A. 
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We  have  looked  at  the  3-wave  instability  for  the  parameters  of  the  UCI 
stellatron  (see  Fig.  14),  and  find  that  growth  is  predicted  in  the  energy 
range  2. 5-3. 5  MeV,  with  the  intersection  in  Fig.  15  occurring  at  about  t  = 
18.  The  predicted  growth  rate  is  quite  large,  about  30  e-foldings/ usee. 

At  this  point  we  do  not  know  whether  any  of  the  current  loss  in  the  experi¬ 
ment  can  be  attributed  to  this  instability.  Thermal  effects  have  not  been 
included  in  the  theory,  and  at  high  t-numbers  may  stabilize  the  mode.  A 
quadrupole  focusing  experiment  being  mounted  at  PSI  has  been  deliberately 
designed  to  push  the  intersection  in  Fig.  15  out  to  high  t-numbers  [11]. 


IV.  SOLENOIDAL  LENS  BETATRON 


1.  BEAM  TRANSPORT 

The  solenoidal  lens  betatron  is  a  strong-focusing  device  with  sole¬ 
noids  of  alternating  polarity  placed  around  a  torus  [4].  As  the  beam 
travels  around,  it  rotates  poloidally  at  the  Larmor  frequency,  changing  its 
direction  of  rotation  from  one  solenoid  to  the  next.  A  device  of  this  type 
has  been  constructed  at  the  Institute  for  Accelerator  and  Plasma  Beam 
Technology  (IAPBT)  at  the  University  of  New  Mexico. 

We  have  simulated  beam  transport  in  a  solenoidal  lens  betatron  using 
both  the  exact  fields  for  finite-length  solenoids,  and  a  sinusoidal 
approximation  to  the  fields,  with  similar  results.  The  simulations, 
described  in  Appendix  8,  demonstrate  the  higher  mismatch  tolerance  of  the 
device  relative  to  the  modified  betatron.  Recent  experiments  on  the 
machine  have  convincingly  shown  this  property  by  being  able  to  transport 
current  around  180  degrees  without  any  vertical  field  [4],  This  tolerance 
is  considerably  larger  than  that  seen  in  the  simulations  in  Appendix  B 
because  the  experimental  injection  energy  is  just  300  kV  instead  of  1  MeV, 
so  that  that  the  centrifugal  force  on  the  beam  is  substantially  smaller. 

If  this  large  energy  acceptance  is  to  be  maintained  as  the  beam  is 
accelerated,  then  the  solenoidal  fields  must  be  ramped  In  synchronism  with 
the  beam  energy.  This  can  be  seen  from  the  expression  for  the  momentum 
compaction  factor 

(7) 

where  Bq  Is  the  amplitude  of  the  solenoidal  field  and  Bz  is  the  vertical 
field.  This  equation  can  be  simplified  if  we  assume  that  the  beam  radius 
Is  matched  to  the  solenoidal  field,  and  that  the  beam  is  space-charge 
dominated,  l.e.,  the  emittance  is  negligible  in  the  force  balance  equation. 
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Since  the  beam  is  injected  from  a  cathode  shielded  from  magnetic  fields, 
then  it  is  always  in  the  Brillouin  limit  (assuming  it  is  matched),  so  that 
the  force-balance  equation  is 


1 

1  _ 
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Substituting  this  into  Eq.  (7)  we  obtain 


(8) 


(9) 


Since  Bz  is  proportional  to  the  beam  energy,  the  compaction  factor  rapidly 
approaches  the  weak-focusing  limit  if  Bg  is  not  increased.  In  Eqs.  (7)-(9) 
we  have  assumed  1  -  n  =  1/2  where  n  is  the  external  field  index.  However, 
the  the  external  index  in  the  IAPBT  betatron  is  essentially  equal  to  unity, 
as  will  be  shown  below.  Thus,  in  the  weak  focusing  limit,  there  is  no 
transverse  focusing  on  the  beam.  It  will  therefore  be  necessary  to  provide 
some  significant  solenoidal  focusing  throughout  the  acceleration.  One 
important  benefit  of  keeping  the  ratio  Bg/Bz  constant  Is  that  no  orbital 
resonances  are  crossed.  Such  resonances  were  a  severe  problem  in  previous 
solenoidal  lens  betatron  experiments  by  dePackh  [12], 


We  can  show  geometrically  that  the  external  field  index  in  the  IAPBT 
betatron  is  effectively  unity.  The  vertical  field  coils  for  the  experiment 
are  wound  on  cylindrical  mandrels  as  shown  in  Fig.  16.  Unlike  the  usual 
sector  magnet  [13],  the  ends  are  not  canted.  From  Fig.  16,  we  see  that  the 
requirement  for  an  equilibrium  orbit  is 
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where  91n,  0out  are  the  entering  and  exiting  angles  of  the  orbit  at  each 
magnet,  and  8S  Is  the  angle  from  the  center  of  one  magnet  to  the  center  of 
the  next.  This  means  that  the  radius  of  curvature  of  a  particle  inside  the 
magnets  must  be 
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where  d  is  the  length  of  the  magnets.  The  vertical  magnetic  field  has  been 
experimentally  measured  to  be  uniform  over  a  radial  interval  of  ±2  cm  from 
the  axis  [13],  Thus,  for  a  given  particle  energy,  the  equilibrium  orbit 
conditions  (10)  and  (11)  can  be  satisfied  for  a  continuum  of  radial 
particle  positions.  Intuitively,  this  follows  from  the  fact  that  the  path 
length  inside  the  magnetic  field  is  the  same  for  all  orbits  satisfying 
Eq.  (10),  regardless  of  the  length  of  the  orbit  circumference  (see 
Fig.  16). 

2.  BEAM  STABILITY  IN  THE  IAPBT  BETATRON 

A  general  treatment  of  stability  of  the  solenoidal  lens  betatron  is 
complicated  by  the  lack  of  toroidal  symmetry  in  the  equilibrium.  One  can 
overcome  this  problem  in  dealing  with  toroidal  mode-numbers  lower  than  that 
associated  with  the  solenoidal  fields  by  an  averaging  procedure.  Details 
of  this  calculation  are  given  in  Appendix  C.  One  finds  that  the  net  effect 
of  the  solenoidal  lenses  is  to  increase  the  radial  and  vertical  focusing  on 
the  beam  compared  to  a  conventional  betatron.  Qualitatively,  however,  the 
beam  behavior  is  similar  to  the  conventional  betatron  in  that  the  trans¬ 
verse  oscillations  are  uncoupled.  In  the  negative-mass  simulations 
described  below,  we  observe  little  motion  along  the  z  direction,  unlike  the 
modified  betatron. 

The  analytic  model  in  Appendix  C  predicts  a  finite  negative-mass 
transition  energy  due  to  the  strong-focusing,  as  in  the  stellatron  [9], 

Also  like  the  stellatron,  the  strong  focusing  suppresses  the  peak  growth 
rate  of  the  Instability.  This  Is  Illustrated  In  Fig.  17,  which  shows  the 
t  *  5  growth  rates  for  parameters  close  to  those  of  the  IAPBT  betatron. 
These  curves  were  obtained  from  the  dispersion  relation  of  Ref.  1  by  set¬ 
ting  Bft  *  0,  and  by  adding  the  smoothed  focusing  effect  of  the  solenoidal 

0  2 
lenses  into  the  expression  for  uy: 
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Figure  17.  Growth  r*tes  of  the  l  =  5  negative  mass  instability  in  the  UNM 
solenoidal-lens  betatron. 
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(The  vertical  focusing  also  increases,  but  the  vertical  oscillations 
decouple  for  B.  =  0.)  Note  that  if  the  solenoidal  focusing  strength  is 

o 

kept  at  its  initial  high  value  by  ramping  the  solenoidal  field  with  the 
vertical  field,  then  the  peak  growth  rate  is  reduced  by  about  a  factor  of 
five  from  the  case  where  the  solenoidal  field  is  held  fixed. 


We  have  run  a  number  of  simulations  to  compare  with  the  analytic 
predictions.  A  sinusoidal  approximation  to  the  reversing  toroidal  field 
was  used,  with  the  wavelength  chosen  to  correspond  to  30  solenoids  around  a 
1  m  radius  torus.  In  the  IAPBT  betatron  there  are  40  solenoids,  but  10  of 
these  are  on  straight  sections  which  we  do  not  simulate  here.  First,  we 
looked  at  the  l  -  5  mode  on  a  1  MeV,  100  A  beam.  The  beam  radius  was 
1.75  cm  and  the  pipe  radius  was  4.5  cm.  The  peak  solenoid  field 
(134  gauss)  was  chosen  to  have  a  matched  beam  at  1  MeV.  Essentially  no 
growth  of  the  instability  was  seen  in  a  period  of  230  ns.  The  theoreti¬ 
cally  predicted  growth  rate  is  7.6  x  10"^  cm'1,  which  would  produce  about 
5.3  e-foldings  in  230  ns.  (This  growth  rate  is  larger  than  those  in 
Fig.  17  because  of  the  lower  solenoidal  field.)  We  attribute  the  observed 
stability  to  the  spread  in  particle  circulation  frequencies.  Since  the 
beam  has  no  net  poloidal  rotation  due  to  the  alternating  direction  of  the 
toroidal  magnetic,  Eq.  (3)  is  always  satisfied.  Thus  Eq.  (2)  should  be  an 
approximate  stability  criterion,  as  it  is  for  the  modified  betatron  (see 
Sec.  II).  Applying  Eq.  (2)  to  the  present  simulations,  we  obtain  8.75  x 
10'4  cm'1  for  the  right-hand  side,  which  is  marginally  sufficient  for 
stability. 

Assuming  this  is  the  correct  reason  for  the  lack  of  growth  in  the 
simulation,  it  is  clear  that  finite  beam  radius  is  an  important  stabilizing 
mechanism  in  the  IAPBT  betatron.  However,  the  beam  radius  shrinks  as  the 
beam  is  accelerated,  particularly  if  the  solenoidal  field  is  ramped  with 
Bz»  For  the  latter  case,  the  beam  radius  is  approximately 
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Therefore,  we  must  trade  off  the  increased  stability  from  the  strong 
focusing  with  the  decreased  stability  due  to  shrinkage  of  the  beam  radius. 
For  the  case  in  Fig.  17,  for  example,  Eq.  (13)  predicts  that  the  beam 
radius  will  shrink  from  2  cm  at  y  =  1.6  to  0.045  cm  at  y  -  20.  From 
Eq.  (2),  we  find  that  the  t  -  5  mode  is  still  just  marginally  stable  at 
Y  =  20. 

The  smooth  approximation  to  the  focusing  force  becomes  dubious  as  the 
mode-number  l  approaches  Ns / 2 ,  where  N$  is  the  number  of  solenoids.  How¬ 
ever,  since  no  approximation  is  made  in  the  simulations,  such  modes  can  be 
looked  at  numerically.  For  the  same  beam  parameters  as  above,  we  simulated 

the  l  =  15  (=  N  /2)  mode,  and  observed  a  strong  instability,  with  a  growth 
-3  -1 

rate  of  2.2  x  10  cm  .  A  particle  plot  from  the  nonlinear  stage  of  the 
instability  (Fig.  18)  shows  standard  negative  mass  behavior,  i.e.,  radial 
kinking  accompanied  by  toroidal  clumping.  The  instability  has  almost  satu¬ 
rated  at  the  point  reached  in  Fig.  18,  but  1/6  of  the  current  has  been 
scraped  off  onto  the  wall.  If  the  wall  radius  had  been  larger,  the  insta¬ 
bility  might  have  saturated  without  current  loss. 

If  we  naively  apply  the  analytic  dispersion  relation  to  the  l  =  15 

-3  -1 

mode,  we  obtain  a  growth  rate  of  2.2  x  10  cm  ,  the  same  as  the  numerical 

.3 

growth  rate.  However,  the  right  hand  side  of  Eq.  (2)  gives  2.6  x  10 
cm-1,  thereby  predicting  stability  due  to  thermal  effects.  Thus,  it 
appears  that  the  cold  beam  growth  rate  for  1  =  15  is  significantly  larger 
than  that  predicted  by  the  analytic  theory. 

With  respect  to  the  negative-mass  transition  energy,  on  the  other 
hand,  the  prediction  of  the  analytic  theory  does  seem  to  hold  even  at 
l  =  15.  We  have  simulated  this  mode  at  y  =  1.6,  which  is  the  injection 
energy  for  the  IAPBT  betatron,  and  observe  no  instability.  The  analytic 
theory  predicts  stability  at  least  up  to  y  =  2.2  (cf.  Fig.  17). 
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Figure  18.  Particle  plot  taken  during  the  nonlinear  stage  of  the  t 
negative-mass  instability  for  UNM  betatron  parameters. 


3.  HIGH  CURRENT  BEAM  STABILITY 


If  the  IAPBT  experiments  at  100  A  are  successful,  then  higher  current 
experiments  will  be  undertaken.  We  have  performed  a  number  of  simulations 
at  10  kA  to  investigate  negative-mass  instability  growth  in  this  regime. 

In  each  of  the  simulations,  we  set  the  beam  radius  to  2  cm  in  order  to  have 
reasonably  good  resolution  of  the  beam  on  a  50  x  50  mesh.  This  means  that 
the  solenoidal  magnetic  field  necessary  to  get  a  matched  beam  decreases  as 
the  energy  increases.  Thus  the  linear  growth  rates,  shown  in  Fig.  19,  do 
not  represent  those  experienced  by  an  accelerating  beam  where  the  solenoi¬ 
dal  field  is  kept  constant,  or  ramped  up.  We  assumed  20  solenoids  around 
the  torus,  and  simulated  the  mode  l  =  N$  =  20.  In  Fig.  19,  we  compare  the 
simulation  results  with  (a)  growth  rates  obtained  for  a  conventional  beta¬ 
tron  with  B0  =  0,  and  (b)  growth  rates  obtained  by  applying  the  smoothed 
approximation  for  the  solenoidal  lenses.  As  in  the  100  A  simulations 
described  above,  we  find  that  the  high-*  growth  rates  are  anomalously  large 
when  thermal  effects  are  taken  into  account.  Nevertheless,  the  simulations 
for  y  *  7  and  y  *  12  both  saturate  in  a  manner  similar  to  the  100  A  case 
shown  in  Fig.  18,  with  little  loss  of  current.  Plots  from  the  nonlinear 
stage  of  these  simulations  are  shown  in  Fig.  20.  There  is  considerable 
churning  of  the  beam  after  these  plots  are  taken,  and  this  stage  of  the 
instability  is  probably  not  well  represented  by  the  single-mode  simulation 
scheme.  The  simulation  for  y  3  17  was  not  run  long  enough  to  see 
saturation. 

Experimental  results  from  the  UCI  stellatron  suggest  that  development 
of  the  negative-mass  instability  does  not  destroy  the  beam,  and  that  most 
of  the  current  can  be  accelerated  beyond  the  point  where  the  instability 
saturates.  The  above  simulations  suggest  that  this  may  also  be  true  for 
the  solenoidal  lens  betatron. 
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Figure  19.  Comparison  of  simulation  growth  rates  for  the  I  -  20  mode  on  a 
10  kA  beam  with  analytic  growth  rates  obtained  for  a  conven¬ 
tional  betatron  (dashed  line)  and  by  extrapolation  of  the 
smoothed  focusing  result  (solid  line). 
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Figure  20.  Nonlinear  development  of  the  £  =  20’  negative  mass  instability 
on  a  10  kA  beam  in  a  solenoidal  lens  betatron.  In  (a),  the 
beam  energy  is  y  *  7  and  in  (b),  y  =  12. 
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V.  CONCLUSIONS 


We  have  Investigated  thermal  effects  on  the  negative-mass  instability 
in  the  modified  betatron.  For  a  monoenergetlc  beam,  the  stabilization  is 
strongly  dependent  on  the  poloidal  rotation  frequency  of  the  beam.  When 
this  frequency  is  close  to  zero,  we  observe  total  stabilization  of  the 
instability  at  high  mode-numbers. 

For  parameters  close  to  the  NRL  betatron  at  1  kA,  our  calculations 
agree  with  the  observed  stability  of  the  beam  at  the  injection  energy.  We 
predict  linear  instability  at  low  toroidal  mode-numbers  as  the  beam  is 
accelerated  above  1  MeV.  The  nonlinear  development  of  the  instability  is 
unknown  for  a  1  kA  beam.  Simulations  of  10  kA  beams  show  loss  of  50%  or 
more  of  the  beam  current  at  low  energies.  Several  kiloamperes  remain 
circulating  after  the  instability  has  saturated,  so  that  a  1  kA  beam  may  be 
able  to  survive  without  current  loss,  although  deterioration  of  beam 
quality  is  to  be  expected. 

We  have  simulated  beam  stability  at  high  energies  and  high  t-numbers, 
in  the  regime  of  the  so-called  hybrid  instability.  In  all  such  simulations 
for  10  kA  beams,  we  observe  saturation  of  the  instability  without  loss  of 
beam  current.  There  is  much  less  transverse  kinking  of  the  beam  than  in 
the  negative-mass  instability  regime  which  occurs  at  lower  energies. 

Calculations  for  the  UCI  stellatron  show  that  the  presence  of  back¬ 
ground  ions  has  a  larger  stabilizing  effect  on  the  negative-mass  insta¬ 
bility  than  the  helical  quadrupole  field.  The  predicted  energy  for  onset 
of  the  negative-mass  instability  is  in  agreement  with  experimental  observa¬ 
tions.  Most  of  the  beam  current  is  observed  to  survive  the  disruption  due 
to  the  instability,  but  the  effect  on  beam  quality  has  not  been  measured. 

A  three-wave  electromagnetic  interaction  was  discovered  in  the 
stellatron.  For  high  currents  and  strong  focusing,  its  growth  rate  can  be 
very  large.  However,  it  is  possible  to  choose  parameters  such  that  the 


interaction  is  pushed  out  to  large  mode-numbers,  where  thermal  effects  may 
stabilize  it. 

A  theory  of  the  negative-mass  Instability  in  the  solenoidal  lens 
betatron  was  derived  using  a  smooth  approximation  for  the  focusing  effects 
of  the  solenoidal  lenses.  The  theory  predicts  stability  for  the  100  A  UNM 
betatron  from  its  injection  energy,  y  =  1.6,  to  y  *  2.2.  Above  this 
energy,  simulations  show  a  strong  negative-mass  instability,  which  produces 
a  transverse  kink  on  the  beam  with  an  amplitude  comparable  to  the  wall 
radius.  Some  current  loss  results.  Simulations  with  currents  up  to  10  kA 
show  similar  results,  suggesting  that  most  of  the  current  can  survive  the 
instability. 
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APPENDIX  A 


ELECTROMAGENTIC  INSTABILITY  IN  A  QUADRUPOLE-FOCUSING  ACCELERATOR 


Electromagnetic  instability  in  a  quadrupole-focuslng  accelerator 
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Mission  Research  Corporation.  1 720  Randolph  Road.  S.  £..  Albuquerque.  New  Mexico  87106 
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The  addition  of  helical  quadrupole  focusing  to  a  modified  betatron  configuration  is  shown  to  give 
rise  to  an  electromagnetic  instability  under  certain  conditions.  The  instability  arises  from  a  three- 
wave  coupling  between  the  helical  field,  a  transverse  mode  on  the  beam,  and  a  transverse-electric 
waveguide  mode.  An  analytic  dispersion  relation  is  derived.  Several  features  of  the  instability  are 
confirmed  using  three-dimensional  computer  simulations. 


I.  INTRODUCTION 

Betatrons  and  other  recirculating  accelerator  designs 
have  been  studied  in  recent  years  as  compact  accelerators  for 
high-current  electron  beams.1*6  Conventional7  and  modi¬ 
fied1  betatrons  require  precise  matching  of  the  beam  energy 
to  the  vertical  magnetic  field  if  a  fixed  major  radius  is  to  be 
maintained.  To  overcome  this  restriction,  strong  focusing  in 
the  form  of  a  helical  magnetic  quadrupole  has  been  added  to 
the  modified  betatron  configuration. 3*  This  considerably  in¬ 
creases  the  tolerance  of  the  device  to  mismatch.  In  this  pa¬ 
per,  we  show  that  helical  quadrupole  focusing  can  lead  to  an 
electromagnetic  instability  arising  from  a  three-wave  inter¬ 
action  between  the  static  helical  field,  a  transverse  mode  on 
the  beam,  and  a  transverse-electric  (TE)  waveguide  mode. 
This  behavior  is  reminiscent  of  that  which  occurs  in  planar 
geometry  in  the  presence  of  a  rippled  magnetic  field.’  Unlike 
the  free-electron  laser  instability,  which  was  recently  ana¬ 
lyzed  in  the  presence  of  helical  quadrupole  focusing, 10  longi¬ 
tudinal  bunching  of  the  beam  plays  no  significant  role  in  the 
instability  that  we  describe  here.  Also,  the  instability  is  es¬ 
sentially  independent  of  the  radius  of  curvature  of  the  de¬ 
vice.  It  is  thus  not  related  to  the  negative-mass  instability,  for 
example.  The  main  effect  of  finite  radius  of  curvature  is  to 
discretize  tne  toroidal  mode  numbers.  This  can  be  impor¬ 
tant.  since  if  the  instability  width  is  narrow  enough,  it  can  fall 
between  two  allowed  mode  numbers  and  disappear  (sew  Sec. 
II). 

In  Sec.  II,  we  give  an  analytic  theory  of  the  instability.  In 
Sec.  Ill,  we  obtain  a  simplified  dispersion  relation  that  yields 
an  instability  criterion.  Section  IV  details  comparisons 
between  the  analytic  theory  and  three-dimensional  particle 
simulations.  Section  V  gives  our  conclusions. 

In  writing  down  equations,  we  use  convenient  units, 
where  the  electron  charge  e,  mass  m,  and  velocity  of  light  c 
are  scaled  out.  Lengths  are  normalized  to  c/<u0a*  1  cm,  fre¬ 
quencies  to  wo,  velocities  to  c,  densities  to  wlm/im2,  and 
electric  and  magnetic  fields  to  met 0o/e.  To  get  frequencies 
and  growth  rates  in  sec  ~ 1 ,  multiply  the  normalized  values  by 
3X10'°. 

II.  ANALYTIC  MODEL  OF  INSTABILITY 

The  geometry  and  coordinate  system  we  use  are  illus¬ 
trated  in  Fig.  1.  An  electron  ring  of  major  radius  R  and 
minor  radius  a  circulates  around  the  center  of  a  conducting 
torus  whose  major  and  minor  radii  are  R  and  b,  respectively. 


The  externally  applied  magnetic  fields  consist  of  a  vertical 
field  where  s  is  the  so-called  field  index  ( we  assume 

s  =  1  /2),  a  solenoidal  field  Bt,  and  a  helical  quadrupole  field 
B, .  This  is  the  configuration  of  the  stellatron  accelerator. 3  A 
configuration  where  the  helical  quadrupole  field  is  generated 
by  two  current-carrying  wires  instead  of  four  has  also  been 
proposed.6  We  assume  that  the  beam  can  be  modeled  as  a 
string  of  rigid  disks.  This  means  that  we  follow  the  trans¬ 
verse  motion  of  the  beam  centroid,  and  ignore  any  internal 
dynamics.  This  treatment  is  valid  provided  the  drift-tube 
minor  radius  is  much  greater  than  the  beam  radius,  b>a, 
since  the  fields  we  need  to  consider  have  a  transverse  scale 
length  ~b.  We  ignore  the  perturbed  0  motion  of  the  beam 
under  the  assumption  of  relativistic  stiffness  in  this  direc¬ 
tion,  i.e.,  » 1,  where  y  is  the  beam  relativistic  factor. 

In  equilibrium,  the  matched  value  of  the  vertical  mag¬ 
netic  field  is  B,  —  yVg/R,  where  Ve  is  the  toroidal  beam 
velocity.  Linearizing  about  the  equilibrium  position,  we  ob¬ 
tain  the  following  equations  of  motion  for  the  perturbed  co¬ 
ordinates  r,  z  of  the  beam  centroid3: 

r  +  j  il]f  —  n 9z  +  fstl](r  cos  md  +  z  sin  md) 

=  (-Er-VtB')/y,  (1) 

z  + 1  tl\z  +  tig r  +  fitll  (f  sin  md  —  z  cos  md) 

**{-E,  +  VtBr)/y,  (2) 

where  a  dot  denotes  the  total  time  derivative  d/dt 
-I-  ( Vg/R  )d  /dd,  ft,  =  B,/y,  ft,  =  Be/y,  E  and  B  denote 
the  perturbed  se//-electric  and  self-magnetic  fields  acting  on 


FIO.  I .  Geometry  of  stellatron  accelerator  showrni  conduct m|  toroidal 
cavity  and  externally  applied  magnetic  fields. 
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the  beam  at  the  perturbed  position,  and  p  denotes  the  ampli¬ 
tude  of  the  quadrupole  field  index.  The  quadrupole  field  in¬ 
dex  is  a  normalized  measure  of  the  transverse  gradient  of  the 
quadrupole  field.  It  is  defined  by  (R  /B, )  (dBq,/dr)  evaluat¬ 
ed  at  r  =  R,  and  has  the  form  p  cos  md.  Letting  £  =i  +  i  r, 
Eqs.  ( 1 )  and  (2)  can  be  written  as 

=  -  [E,+iEr+iVe(B,+iBr)]/Y,  (3) 

where  *  denotes  the  complex  conjugate.  To  calculate  E  and 
B,  we  use  a  Green’s  function  approach.  Maxwell’s  equations 
can  be  written  as 

VXVxE-<u2E  =  /«J.  (4) 

B=  -  {i/(o)VxE,  (5) 

where  we  assume  an  e"“  dependence  for  the  perturbed 
fields  and  currents.  The  solution  to  Eq.  (4)  can  be  written 
down  by  using  a  dyadic  Green’s  function,"  constructed 
from  the  solutions  to  the  homogeneous  counterpart  of  Eq. 
(4).  The  expressions  obtained  are  infinite  series."  In  this 
paper,  we  will  deal  only  with  the  interaction  of  the  beam  with 
the  lowest  electromagnetic  waveguide  mode,  i.e.,  the  TE, , 
mode.  Therefore,  we  extract  from  the  infinite  series  for  the 
fields  those  terms  containing  the  TE,,  resonance.  All  other 
self-field  contributions  are  neglected.  This  procedure  is  equi¬ 
valent  to  treating  the  beam  as  a  weak  perturbation  of  the 
vacuum  TE,,  eigenmode.  For  example,  the  expression  for 
the  eM  component  of  £,  evaluated  at  the  location  of  the  beam 
is 


-pnU  £,-)•  =  o, 


Er  =  —  v<d(<u  —  /fl,)? 


v2u 


+  l2/R 2  -  co2 


+  0 


(*)• 

(6) 


where 

l£'  =  {(b2  -Vu'UfiVnb) 

( where  7,  is  a  Bessel  function ) ,  77 , ,  denotes  the  lowest  root  of 
7 ;  ( 7/b )  =  0,  and  0(b/R)  denotes  toroidal  corrections, 
which  we  neglect  under  the  assumption  R>b.  Also,  v  is 
Budker’s  parameter,  i.e.,  the  number  of  electrons  per  unit 
length  of  the  beam  times  the  classical  electron  radius.  In 
deriving  Eq.  (6),  we  have  used  the  rigid-disk  approximation 
to  write  J.  =  —pf  =  i(co  —  /ft, )  pf,  where  p  is  the  beam 
density.  Similar  expressions  are  obtained  for  the  other  fields. 

We  now  proceed  to  solve  Eq.  (3)  by  writing  £  in  the 
form 


'+  gr*"* 


(7) 


where  + ,  —  refer  to  “forward"  and  “backward”  waves, 
respectively.  Substituting  this  into  Eq.  (3),  we  obtain  the 
coupled  equations 


(-n-.+lni-n.n ,  + 


( 


-pn2a  :„_,)•  =  o, 

-  (ft*  -  mil, ):  +  yft,J  +  (ft*  -mft,)fl9 


(8a) 


g(fU  -  md,)2 
co2  -  -  (I  +  m)2/R 


t) 


f-fll  +±n>  +  n_ns+ - — - ) 

\  2  u2-rn  -1-/R-) 

x£r  -pClU  £!„_,)•  =  0, 

(H_  yfl*-  (n-  +mftl)ne 


(8b) 


(8c) 


-I-  • 


a(ft_  +  mft, ) 2 


0 


(9) 


co2  -  7]2n  -  U  +  m )  2/R  2 , 

=  0.  (8d) 

where  fl  ±  =  to  ^  Id,  and  a  =  vlv  /y~  iAv/yb 2.  Note  that 
this  is  a  closed  system  of  equations,  in  which  only  mode 
numbers  /  and  —  m  —  l  appear.  This  is  a  consequence  of  the 
complex-conjugate  sign  appearing  in  Eq.  ( 3 ).  Further,  note 
that  £  *  is  coupled  only  to£  Im_,.  and  £ is  coupled  only 
t0  £  -  m  -  /•  Thus,  the  4  X  4  matrix  of  Eq.  ( 8 )  splits  into  two 
2x2  matrices.  The  determinant  of  one  2x2  matrix  gives  the 
dispersion  relation 

(nl  ~Ttt‘+a*a,~^-n2i-i2/R2) 
xun.  —  mfl.j’-i-  n;  -  <n.  -mn.in, 

- a(n.  \ 

OJ2  -  T)\x  -  (/  +  m)‘/R  2  J 

The  determinant  of  the  other  2x2  matrix  can  be  obtained 
from  Eq.  (9)  by  letting  co-*  —  co.  Having  obtained  a  root  for 
co  from  Eq.  ( 9 ) ,  the  corresponding  normal  mode  for  £  can  be 
obtained  from  Eqs.  (8a)  and  (8b),  and  is  a  linear  com¬ 
bination  of  two  terms  with  space-time  dependencies 
exp(i 16  —  icot)  and  exp[  -  i(m  +  1)9  +  icot  ],  respectively. 

Multiplying  out  the  dispersion  relation  in  Eq.  (9),  we 
obtain  an  eighth-degree  polynomial  in  co.  Numerical  solu¬ 
tion  for  the  sample  parameters  in  Table  I  gives  the  set  of 
curves  shown  in  Fig.  2.  The  parameter  e  ( which  is  frequently 
used  by  stellarator  physicists)  in  Table  I  is  related  to  p  by 
p  =  emilg/ 20,.  The  curves  can  be  classified  as  follows: 

Electromagnetic  modes: 

co2xT)2tl  +  l2/R2,  +  («  +  D2/R 2 , 

Cyclotron  modes: 

cocz(m  +  /)n,  +  cig,  in,  -  ntf . 

Low-frequency  transverse  modes: 

coz:(m  +  l)n,  —  coM,  /fl,  ■+■  cot  , 

TABLE  I.  High -current  stellatron  parameters  used  in  particle  simulations 


(10) 


Torus  major  radius 
Torus  minor  radius 
Beam  major  radius 
Beam  minor  radius 
Beam  current  II) 

Beam  energy  ( y) 

Toroidal  magnetic  Aeld 
Vertical  magnetic  Held 
Quadrupole  Aeld  (e) 
Quadrupole  mode  number  (m) 


I  m 

8. 4-9.5  cm 

1  m 

2  cm 

300  A-IOkA 
7 

5kG 
1 18  G 
0.7 
14-30 
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FIG  2  Frequencies  and  growth  rates,  denoted  by  <u,  and  I\  respectively, 
obtained  from  Eq  (9)  for  Table  I  parameters  with  kA. 

m  a  22  In  (a),  the  two  unstable  interactions  are  circled.  In  (b),  growth 
rates  at  integer  ( allowed )  values  of  /  are  shown.  Multiply  ai„  T  by  3  X  10'° 
to  get  values  ui  units  of  sec" '. 

- - 1 


where  aB  s  l/4ft2/09 .  These  approximate  forms  are  based 
on  the  assumption  ft9  >ft, ,  which  is  satisfied  for  typical  stel- 
latron  parameters. 

It  can  be  shown  from  Eq.  ( 9 ) ,  and  is  apparent  from  Fig. 
2,  that  the  dispersion  relation  is  symmetrical  about  the  line 
I  =  —  m/2  (m  is  an  even  integer,  equal  to  twice  the  number 
of  minor  turns  the  helical  conductors  make  in  going  one 
major  tum  around  the  torus).  More  precisely,  if  we  define 
/'  =  /  +  m/2,  then  a— ►  —  a*  as  /'—►  —  /'.  In  the  following, 
we  will  look  only  at  / '  >  0.  Results  for  / '  <  0  follow  by  sy  mme- 
try. 

Instability  occurs  as  a  result  of  the  intersection  of  the 
curve  as  (m  +  l)fl,  —  ud  with  the  electromagnetic  mode 
as  (irf,  4-  l2/R  2)'n.  The  instability  can  be  thought  of  as  a 
parametric  process'2  in  which  the  quadrupole  field  plays  the 
role  of  a  pump  wave.  If  we  denote  the  three  interacting 
modes  by  subscripts  1,2,3  then  we  get  the  following  frequen¬ 
cies  and  mode  numbers: 

Pump  wave: 

a,  =  0,  k,  =  m/R  , 

Electromagnetic  wave: 

a3  =  (V2,  +12/R  2)'n,  *,-//*,  (11) 

Slow  transverse  wave: 

o>,=  (m  +  /)ft„  k,=  (m +  /)//(, 

where  we  have  neglected  the  small  aB  term  in  a}.  In  terms  of 
this  picture,  instability  occurs  when  the  usual  matching  con¬ 
ditions  are  approximately  satisfied  (a,  4-  a2  = 
k,  +  k2  =  kj).  Energy  for  the  instability  is  provided  by  the 
beam  via  the  negative-energy  slow  transverse  mode. 

III.  DERIVATION  OF  INSTABILITY  CRITERION 

An  approximate  expression  for  the  growth  rate  can  be 
obtained  by  letting  a  =  an+6,  where  wn  =  (77}, 
+  l2/R  2)'/2,  in  Eq.  (9),  and  keeping  terms  to  order  <52.  This 
gives  the  quadratic  expression 


2eW|  1  Cct*  1  ]  —  /ft, )  (d>n  —  /ft,  +  ft#  )ft»52  4-  {2<u  1 1  (<r>  1 1  —  /ft, )  (at  1 1  —  /ft,  4-  ft9 )  [  (  /  +  m)ft,  —  o>,  1  ]  ft9 

+  (<u,,  —  /ft,  )2aft#  -  2a>,,/f2ft;}6  -  a(a»,,  -/ft,)2  [(m-|-/)n,  -<u,,]ft9  =0. 


(12) 


Defining  the  quantities  y, ,yi,y}  by 


y,  **  2a,,(an  -  ia,)(an  -  /a,  +  ft#) 

X[(/  +  m)ft,  -a>n]ft,. 

(13) 

Yi  m  ®(a>n  —  /ft,  )20#  ,  y}  -  2v,,M2(li  , 
the  condition  for  instability  can  be  written  as 

(Vi  +  r2-yi)2<4y,y2, 

(14a) 

or,  equivalently, 

(Vi  -r2-yj)2<4y2y3. 

(14b) 

Equations  ( 12)— ( 14)  allow  us  to  make  some  more  exact 
statements  about  the  conditions  for  instability.  Fran  Eqs. 

( 13)  and  ( 14b),  we  see  that  the  instability  disappears  for 
H  »  0,  as  one  would  expect.  Equation  ( 14a)  shows  that 
(/  +  m)fl,  >a>,,  (i.e.,  y,  >0)  is  required  for  instability  (as* 
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suming  ft9>0;  the  case  ft9  <0  is  discussed  later).  This 
means  that  the  instability  turns  on  when  the  frequency  of  the 
transverse  beam  mode  is  somewhat  above  the  TE, ,  eigenfre- 
quency,  rather  than  exactly  equal  to  it.  To  see  how  the  insta¬ 
bility  turns  off  as  /  increases,  we  rewrite  Eq.  ( 14a)  to  obtain 

~  VftlV^uft;  (15) 

as  the  criterion  for  instability.  The  term  y,  is  sensitive  to  the 
value  of  /  through  the  factor  (/  +  m)ft,  -  a,,.  Ml  in¬ 
creases,  y,  eventually  becomes  large  enough  to  violate  this 
inequality,  and  the  instability  disappears  (cf.  Fig.  2).  Equa¬ 
tion  { 15 )  also  shows  that  as  v— 0  ($0  that  y2-*0)  and O, 
the  unstable  region  in  /  space  becomes  narrower.  If  /  were  a 
continuous  variable,  then  the  instability  would  persist  near 
(/  +  m)ft,  =  au  as  long  as  v  and  n  were  finite.  Since  /  is 
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discrete  in  a  toroidal  system,  however,  the  unstable  region 
can  fall  between  two  integer  values  of  /,  and  no  instability 
would  be  seen. 

Finally,  from  Eqs.  (12)  and  ( 14a ) ,  the  peak  growth  rate 
as  a  function  of  n  is  found  when  /i  is  chosen  so  that 
Y\  +  Yt  —  Yi  =  0.  The  growth  rate  T  is  then  given  by 

r2_  a  [  (m  4-  /)fl,  — »n]  (&>,,  —  ICl, ) 

2  -/n, +n#) 

For  the  parameters  of  Fig.  2,  this  equation  predicts 
T  =  1.06x  10-2  ( 3. 18 X  10*  sec-1  in  cgs  units)  for  /  =  0,  in 
good  agreement  with  the  exact  result  of  1.12X  10-2,  which 
occurs  for  e  =  0.8  ( /is 380). 

Thus  far,  we  have  assumed  Cl,  >  0.  If  ft,  is  negative, 
then  Eq.  ( 14b)  shows  that  the  instability  disappears.  This  is 
to  be  expected  since  the  slow  transverse  wave  then  becomes  a 
fast,  positive-energy  wave.  At  the  same  time,  the  cyclotron 
wave  us  ( m  +  /)fl,  +  0#  becomes  a  slow  wave,  with  the 
potential  for  an  unstable  interaction.  We  have  not  examined 
this  case,  however,  since  previous  calculations1  have  shown 
that  the  parameter  space  in  which  single-particle  orbits  are 
stable  shrinks  greatly  when  the  sign  of  Cl,  is  opposite  to  that 
of  m. 

IV.  COMPARISON  WITH  NUMERICAL  SIMULATIONS 

To  verify  the  above  analytic  calculations,  we  have  per¬ 
formed  three-dimensional  numerical  particle  simulations 
with  the  code  ivoky.11  In  ivory,  field  quantities  are  repre¬ 
sented  on  a  spatial  mesh  in  the  r-z  plane,  and  by  Fourier 
modes  in  the  6  direction.  A  given  field  component  thus  has 
the  form 

F(r,0j,t)  =  £/„  (W)cos  nd  +  g„  (r,z,/)sin  nd 

m 

(17) 

The  number  of  Fourier  modes  kept  in  the  sum  depends  on 
the  problem.  At  a  minimum,  we  must  retain  those  modes 
that  the  linear  theory  indicates  are  coupled  together.  In  addi¬ 
tion.  the  n=Q,m  modes  must  be  kept  since  the  equilibrium 
fields  contain  components  in  these  modes.  The  self-consis¬ 
tent  fields  of  the  beam  are  advanced  in  time  using  the  full 
Maxwell  equations.  Particles  are  advanced  using  the  full 
Lorentz  force  equations.  The  stellatron  fields  are  computed 
from  analytic  expressions.1 

In  deciding  cm  parameters  for  the  simulations,  comput¬ 
ing  costs  constrained  us  to  choose  cases  that  minimize  the 
running  time  and  storage  requirements.  Thus,  we  concen¬ 
trated  on  cases  with  large  expected  growth  rates.  In  addition, 
we  either  chose  /  =  0,  or  chose  m  +  /  to  be  a  small  integer 
multiple  of  /.  (Recall  that  /  denotes  the  TE,,  wave  mode 
number. )  This  minimizes  the  number  of  particles  needed  to 
resolve  the  different  mode  numbers.  For  NO,  for  example, 
we  see  from  Eqs.  (11)  and  ( 17)  that  only  mode  numbers  0, 
m  must  be  represented  in  the  simulations. 

Choosing  /  -  0,  we  used  the  parameters  in  Table  I  with 
b  —  8.8  and  a  10  ItA  beam  current,  and  performed  simula¬ 
tions  for  different  values  of  the  quadrupok  mode  number  m. 
The  growth  rates  obtained  are  plotted  in  Fig.  3  versus  those 
obtained  from  Eq.  (9).  We  also  show  the  growth  rate  from  a 
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FIG.  3.  Growth  rates  of  instability  versus  quadrupok  mode  number  m 
(which  must  be  an  even  integer)  for  parameters  m  Table  I.  and  6  —  8  8. 
/  —  10 kA  Theoretical  results  (solid  line)  are  compared  to  simulation  code 
results. 

simulation  for  m  —  22  in  which  the  beam  was  represented  by 
a  string  of  rigid  disks  in  the  code,  instead  of  the  more  realistic 
particle  representation.  There  are  at  least  two  possible  rea¬ 
sons  why  agreement  with  the  analytic  result  is  better  for  the 
rigid  disk  simulation.  First,  the  initial  field  energy  level  is 
lower  for  the  latter  as  a  result  of  the  absence  of  internal  de¬ 
grees  of  freedom  of  the  beam,  so  that  there  is  a  longer  visible 
period  of  linear  growth.  In  the  simulations  with  panicles,  the 
field  energy  increases  only  about  one  order  of  magnitude 
over  its  initial  value  before  saturating,  as  shown  in  Fig.  4. 


TIMS  (wawc) 

FIG.  4.  Energy  in  e  •  22  Bek)  components  ( r.  -  toroidal  mode  number ) 
versus  time  showing  development  of/-  0,  m  »  2,  instability  in  Fig.  2.  This 
plot  was  generated  by  the  simulation  code  ivory. 
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FIG  5.  Contour  plot  of  n  =  0  component  of  toroidal  self-magnetic  field  B, 
during  early  nonlinear  growth  of  an  /  =  0.  m  =  20  instability.  The  contours 
have  linearly  increasing  values  from  A  to  G.  The  structure  is  that  of  a  TEn 
mode  The  distortion  of  contours  in  the  center  is  a  result  of  the  presence  of 
the  beam. 

This  means  that  the  growth  measured  may  not  be  truly  expo¬ 
nential.  Second,  in  the  particle  simulations  there  is  a  spread 
in  the  energy  of  the  particles  as  a  result  of  space  charge, 
giving  rise  to  a  spread  in  transverse  oscillation  frequencies. 
This  may  have  a  stabilizing  effect. 

There  is  good  agreement  between  simulations  and  the¬ 
ory  as  to  the  turn-on  and  tum-off  of  the  instability.  The  TE, , 
cutoff  frequency  for  the  drift  tube  is  77,  ,  ss  1 .84/6  ~  0.2 1 .  Fig¬ 
ure  3  shows  that,  in  agreement  with  the  analytic  prediction 
(Sec.  Ill),  instability  sets  in  for  m  >  77,,/fl,~21.  No  insta¬ 
bility  is  seen  during  the  length  of  the  simulations  for  m  >  26 
(we  ran  cases  for  m  —  28,30). 

Further  evidence  of  the  electromagnetic  character  of  the 
instability  comes  from  two  /  =  0  simulations  where  the  mi¬ 
nor  cross  section  of  the  drift  tube  was  varied,  keeping  m 
fixed.  For  m  =  20,  /  =  10  kA,  6  =  8.8  (other  parameters  as 
in  Table  I),  no  growth  was  observed,  since  mfl,  < tju.  On 
increasing  the  minor  radius  to  9. 5, 7 , ,  decreased  to  give  mfl, 
-i)M,  and  strong  growth  was  observed.  In  addition,  a  con¬ 
tour  plot  of  the  magnitude  of  the  n  =  0  component  of  the 
perturbed  Be,  shown  in  Fig.  5,  shows  the  sin(  ^  4-  d0) 
x/,  (Vt  1 P 1  dependence  of  a  TE, ,  mode.  Here.p,  <t>  are  polar 
coordinates  in  the  plane  of  Fig.  3,  and  ^0is>  phase  factor. 

In  order  to  see  whether  the  aspect  ratio  of  the  torus 
affects  the  instability,  we  increased  the  major  radius  R  first 
to  103  and  then  to  10\  keeping  m/R  fixed  at  0.22  and  choos¬ 
ing  /  =  0.  No  significant  change  in  the  growth  rate  was  ob¬ 
served  in  the  simulations.  This  is  in  contrast  to  the  negative- 
mass  instability,'4  where  the  growth  rate  goes  to  zero  as 
R— »oo  • 

Next,  we  tested  the  vl/2  dependence  of  the  growth  rate 
on  the  beam  current  predicted  by  Eq.  ( 16).  The  parameters 
in  Table  I  with  b  =  9.0,  m  =  22  were  used  to  perform  simu¬ 
lations  at  300  A,  1  kA,  and  10  kA  beam  currents.  The  results 
in  Fig.  6  show  that  the  predicted  scaling  is  supported  by  the 
simulations. 
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FIG.  6.  Growth  rite  of  instability  versus  beam  current  for  Table  I  param¬ 
eters  with  m  =  22,  b  =  9.0.  Simulation  results  are  compared  to  theoretical 
prediction. 


Thus  far,  we  have  reported  on  results  for  1  =  0.  Now,  we 
look  at  a  case  for  which  /  =  m  =  14  The  dispersion  diagram 
for  this  case,  obtained  from  Eq.  (9),  is  similar  to  that  in  Fig. 
2,  except  that  one  of  the  unstable  regions  occurs  around 
/  =  14.  The  simulation  parameters  used  are  those  in  Table  I, 
with  7=10  kA,  b  =  8.4  chosen  to  enhance  the  growth  rate. 
From  Eqs.  (11)  and  ( 17),  we  see  that  mode  numbers  0,  14, 
28  must  be  represented  in  the  simulation  code.  Field  energy 
in  the  modes  n  =  14,  28  is  observed  to  grow  exponentially  at 
a  rate  T  =  6.3x  10“\  compared  to  the  analytic  rate 
7.8  X  10-3.  In  this  case,  it  is  the  n  =  14  field  plots  that  show  a 
TE  mode  character  similar  to  that  in  Fig.  5. 

The  nonlinear  development  of  the  instability  is  an  im¬ 
portant  issue.  To  address  it  completely,  many  modes  would 
have  to  be  kept  in  the  simulation  code,  since  nonlinear  effects 
give  rise  to  the  generation  of  modes  other  than  those  in¬ 
volved  in  the  linear  growth  stage.  However,  we  believe  that 
the  simulation  results  with  just  the  linear  modes  present  may 
give  a  good  guide  to  the  nonlinear  development  of  the  insta¬ 
bility  for  the  following  reasons.  First,  the  quadrupole  field 
gives  an  initial  perturbation  to  the  mode  number  m,  so  that 
the  fields  in  this  mode  are  not  growing  from  random  noise 
(see  Fig.  4).  Second,  although  wave-wave  interactions  are 
not  treated  correctly  if  we  do  not  include  other  modes,  the 
wave-particle  interactions  are  treated  nonlinearly.  The  sim¬ 
ulations  show  that  in  those  cases  with  large  linear  growth 
rates,  which  we  can  therefore  afford  to  run  to  “saturation,” 
the  wave  amplitude  grows  until  the  beam  strikes  the  wall 
(see  Fig.  4) .  As  a  result,  most  of  the  beam  particles  are  lost, 
leaving  a  large-amplitude  TE,,  wave  in  the  drift  tube. 

Finally,  we  look  at  what  our  analysis  predicts  for  the 
stellatron  experiment  presently  under  way  at  the  University 
of  California  at  Irvine  (UCI).1516  The  parameters  for  this 
experiment  are  tabulated  in  Table  II.  Because  of  the  low 
current,  the  width  of  the  instability  is  quite  narrow,  extend¬ 
ing  over  at  most  one  /  number.  As  a  result,  small  changes  in 
parameters  can  cause  the  instability  to  abruptly  appear  and 
disappear  ( cf.  Sec.  Ill ) .  We  find  that  the  beam  is  unstable  in 
the  region  y  =  2.4-2.45.  With  a  one-turn  accelerating  vol¬ 
tage15  of  500  V,  the  beam  would  spend  sufficient  time  in  this 
region  to  undergo  5-6  e  foldings.  This  instability  may  be 
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TABLE  II.  Parameters  of  UCI  stellatron. 


Torus  major  radius 

41  cm 

Torus  minor  radius 

4  cm 

Beam  current 

200  A 

Beam  energy 

20  kV-4  MeV 

Toroidal  magnetic  field 

0-10  kG 

Vertical  magnetic  field 

0-400  G 

Quadrupole  field  (f ) 

0.18 

Quadru  pole  mode-number 

12 

related  to  the  current  disruption  seen  in  some  firings  during 
the  early  part  of  the  acceleration. 14  We  emphasize  that  this  is 
a  tentative  explanation.  It  has  also  been  suggested  that  the 
negative-mass  instability  may  be  responsible  for  this  disrup¬ 
tion.14 

V.  SUMMARY  AND  CONCLUSIONS 

We  have  derived  a  dispersion  relation  for  a  parametric 
electromagnetic  instability  in  a  stellatron  accelerator.  The 
instability  arises  from  the  interaction  between  the  quadra- 
pole  winding,  a  negative-energy  transverse  wave  on  the 
beam,  and  an  electromagnetic  waveguide  mode.  The  growth 
rate  of  the  instability  is  independent  of  the  radius  of  the  to¬ 
roidal  drift  tube.  It  therefore  occurs  in  straight,  as  well  as 
toroidal,  systems.  Three-dimensional  numerical  simulations 
of  the  stellatron  have  been  carried  out  with  the  code  ivory. 
The  simulated  linear  growth  rates  and  conditions  for  the 
onset  of  the  instability  are  in  reasonably  good  agreement 
with  the  analytic  model.  The  simulations  show  strong  dis¬ 
ruption  of  the  beam  in  the  nonlinear  regime,  leading  to  loss 
of  current.  Our  calculations  predict  some  growth  of  the  in¬ 
stability  for  the  parameters  of  the  UCI  stellatron  experi¬ 
ment.  The  instability  may  be  a  factor  in  limiting  the  beam 
current  in  this  experiment. 
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The  solenoldal  lens  betatron1  uses  a  series  of  solen- 
oldal  lenses  arranged  around  a  race-track  shaped 
drift-tube  to  provide  strong  transverse  focusing  for 
a  high-current  electron  ring.  Equilibrium  behavior 
of  the  circulating  beam  Is  examined  for  parameters 
close  to  those  of  the  University  of  New  Mexico  mach¬ 
ine  currently  under  construction.  The  tolerance  of 
the  beam  to  mismatches  in  the  toroidal  and  vertical 
fields  is  evaluated  analytically  and  using  a  particle 
simulation  code.  The  linear  and  nonlinear  develop¬ 
ment  of  the  negative-mass  Instability  In  the  device 
Is  also  studied.  Stability  behavior  comparable  to 
that  in  a  conventional  betatron  (l.e.,  one  with  no 
toroidal  magnetic  field)  Is  found.  Growth  rates  are 
compared  to  those  obtained  from  an  analytic  model. 

Equilibrium  Properties 

In  order  to  confine  and  accelerate  high  currents 
In  a  betatron,  the  conventional  weak  focusing  must  be 
supplemented.  In  the  solenoldal  lens  betatron  ( SLB ) 1 
this  is  accomplished  through  the  use  of  periodic  sol- 
enoldal  lenses,  as  shown  In  Fig.  1.  The  SLB  thus 
differs  from  the  "modified"  betatron,2  for  example, 
which  uses  a  uniform  toroidal  field.  For  a  matched 
equilibrium,  the  beam  In  the  SLB  reverses  its 
pololdal  rotation  at  each  lens,  so  that  the  net 
poloidal  rotation  Is  zero.  This  requires  that  the 
beam  be  Injected  from  a  cathode  which  Is  shielded 
from  magnetic  flux,  so  that  the  beam  produced  has  no 
canonical  pololdal  angular  momentum.  In  the  region 
between  lenses,  the  beam  rotates  at  the  Larmor  fre¬ 
quency,  l.e.,  ne/2v,  where  ne  1s  the  nonrelatlv- 
Istlc  cyclotron  frequency  In  the  solenoldal  field, 
and  y  is  the  relativistic  factor.  If  the  beam  emlt- 
tance  Is  negligible,  then  force-balance  requires 

7  -  ns  +  p2/4  •  0  (1) 

where  ns  *  2vR2/y3B2a2,  Is  the  self-field  index,  v 
Is  Budker's  parameter  (beam  current  divided  by  17 

kA),  R  and  a  are  major  and  minor  beam  radii,  respect¬ 
ively,  e  Is  the  beam  velocity  normalized  to  c,  p  is 
the  ratio  Bg/Bz,  where  Bz  Is  the  vertical 
betatron  field,  and  Bg  Is  the  solenoldal  (toroidal) 
field  between  the  lenses.  We  have  assumed  that  the 
vertical  field  Index  Is  1/2.  From  Eq.  (1),  we  can 
compute  the  solenoldal  field  needed  for  a  matched 

beam,  given  the  other  quantities.  Sample  parameters 
for  Initial  experiments  at  the  University  of  New 
Mexico  (UNM)  are:  R  *  1  m,  a  *  2  cm,  y  «  3,  v  » 
5.9*  10" 1  (100  Amps).  From  these,  we  compute  that 
Bg  ■  80  Gauss  Is  required.  This  result  Is  expected 
to  be  accurate  In  the  limit  vrfiere  the  thickness  of 
the  magnetic  cusps  Is  much  less  than  the  distance 
between  them.  A  scenario  closer  to  that  envisioned 
for  the  UNM  device  Is  to  have  21  cm  long  solenoids 

spaced  10.5  cm  apart.  There  are  thus  20  solenoids 

evenly  spaced  around  the  1  m  major  radius  torus.  We 
model  the  magnetic  fields  by  using  the  exact,  numer¬ 
ically  computed  values  for  finite  length  straight 
solenoids.  One  "cell"  of  the  periodic  field  Is  shown 
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in  Fig.  2.  Simulations  of  beam  behavior  in  this 
device  are  carried  out  using  the  electromagnetic  par¬ 
ticle  code  IVORY.3  We  find  that  to  minimize  the 
envelope  fluctuations,  the  field  In  t(«>_  middle  of 
each  solenoid  must  be  about  100  Gauss.  The  discrep¬ 
ancy  between  this  and  the  analytic  result  Is  presum¬ 
ably  due  to  the  thickness  of  the  cusps,  which  lowers 
the  average  field  value. 


ID) 


ID) 


Figure  1.  Section  of  Solenoidal  Lens  Betatron.  The 
poloidal  velocity  Vpoi  reverses  from  one 
solenoid  to  the  next. 
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Figure  2.  Vector  plot  of  one-half  period  of  toroidal 
field,  showing  location  of  solenoidal 
windings. 


Mismatch  In  Solenoldal  Field 

To  check  the  tolerance  of  the  beam  to  mismatch 
In  the  solenoidal  field,  we  run  two  cases  in  which 
the  value  of  Bg  Is  doubled  and  halved  respectively. 

In  each  case,  we  ensure  that  the  beam  has  zero  canon¬ 
ical  pololdal  angular  momentum.  The  resulting  envel¬ 
ope  oscillation  amplitudes  are  shown  In  Fig.  3.  The 
beam  oscillates  between  radii  of  2  cm  and  0.5  cm  for 
Bg  ■  200  Gauss,  and  between  2  cm  and  3  cm  for  B9 
■  50  Gauss.  Experimentally,  it  should  not  prove  dif¬ 
ficult  to  avoid  these  large  mismatches. 


106. 


(2) 


i(ea) 

Figure  3.  Amplitude  of  envelope  oscillations  for  (1) 
Bg  *  200  Gauss  (inner  bullet),  (11)  Bg 
*  50  Gauss  (dashed  outer  line). 


Figure  4.  Effect  of  1.5*  vertical  field  mismatch  in 
(a)  SL8,  (b)  modified  betatron. 

Mismatch  In  Vertical  Field 

For  a  given  value  of  the  vertical  field  Bz,  the 
matched  major  radius  of  the  beam  Is  r  »  BYC/nz, 
where  r  measures  the  distance  from  the  major  axis  of 
the  torus.  Transverse  oscillations  about  this  radius 
are  determined  from  the  equation. 


•f  *  -  (*- ",  A*)  •*/.*  ■  ^ 

where  b  Is  the  minor  radius  of  the  drift-tube,  and 
4Bz1s  the  mismatch  In  the  vertical  field.  This 
equation  Is  valid  In  the  thln-lens  limit.  In  which 
the  solenoldal  lenses  act  as  a  continuous  radial 
focusing  force  on  the  beam.  Equation  (?)  predicts  a 
momentum  compaction  factor  a  *  (1/2  -  ns  a2/ b2  + 
p2/4)_1.  The  effect  of  a  1.5*  vertical  field  mis¬ 
match  Is  shown  In  the  particle  simulation  In  Fig.  4, 
In  which  we  assume  the  sample  'JNM  beam  parameters 
given  above,  and  Bg  *  100  Gauss.  For  these  para¬ 
meters,  a  «  1,  so  that  the  beam  oscillates  about  an 
equilibrium  major  radius  of  98.5  cm.  For  comparison, 
the  momentum  compaction  factor  of  the  weak-focusing 
modified  betatron  Is  a  *  (1/2  -  ns  a2/b2)-1  *  2.5. 
In  this  case  a  1.5*  mismatch  In  the  vertical  field 
gives  the  beam  an  average  major  radius  of  96  cm.  As 
seen  In  Fig.  4,  the  oscillation  about  this  position 
brings  the  beam  in  contact  with  the  wall.  Even  with 
the  Improved  mismatch  behavior  of  the  SLB,  however, 
keeping  the  vertical  field  matched  to  the  beam  energy 
to  within  a  few  percent  will  require  delicate  tuning 
of  the  experiment. 

Stability  Properties 

A  circulating  high-current  electron  ring  may  be 
subject  to  several  types  of  instabilities,  including 
negative-mass,  resistive  wall,  and.  In  the  case  of 
the  SLB,  which  has  accelerating  gaps,  the  beam  break¬ 
up  Instability.  In  addition,  single  particle  orbital 
resonances  may  affect  beam  quality.  Here,  we  concen¬ 
trate  on  the  negative-mass  Instability,  since  it  is 
potentially  the  most  serious  collective  instability.3 
An  analytic,  high-current  theory  of  this  and  some 
closely  related  Instabilities  Is  described  in  a  com¬ 
panion  paper.4  The  theory  Is  directly  applicable 
only  to  devices  In  which  the  toroidal  coordinate  is 
Ignorable  In  the  equilibrium,  such  as  the  modified 
betatron.  However,  we  use  the  theory  here  as  a  guide 
In  discussing  the  stability  of  the  SLB.  Our  numeri¬ 
cal  results  obtained  from  3-D  simulations  using  IVORY 
are  not  restricted  this  respect.  However,  the  number 
of  simulations  we  can  perform,  and  their  length,  is 
small  due  to  computing  expenses.  We  have  therefore 
concentrated  on  making  runs  at  high  currents,  where 
relatively  large  growth  rates  are  expected. 

The  first  case  we  look  at  is  a  10  kA,  5.5  MeV 
beam  with  a  *  2  cm,  and  R  *  1  m.  We  find  that  a  sol¬ 
enoldal  field  of  600  Gauss  can  transport  this  beam 
with  minimal  ( < 1 0* )  envelope  modulation.  The  solen- 
oidal  lenses  have  the  same  configuration  as  in  Fig. 
2.  The  beam  Is  given  a  small  Initial  perturbation, 
and  the  growth  of  t  *  20  fields  on  the  beam  is  fol¬ 
lowed.  The  justifications  for  following  just  one 
mode  in  the  simulations  are:  a)  In  the  linear  regime, 
mode-coupling  should  be  negligible;  b)  By  following 
one  mode,  direct  comparison  with  linear  theory  pre¬ 
dictions  for  <hat  particular  mode  are  possible;  c) 
Since  there  are  20  solenoids  around  the  torus,  the  t 
»  20  mode  Is  likely  to  have  a  large  Initial  perturba¬ 
tion;  d)  Computing  costs  severely  limit  the  number  of 
toroidal  modes  one  can  afford  to  keep.  In  the  simu¬ 
lation  just  mentioned,  we  find  a  rapid  linear  growth 
rate  of  r  «  1.62*10®  sec*1.  The  effect  of  the  insta¬ 
bility  on  the  beam  as  It  reaches  nonlinear  levels  is 
shown  In  Fig.  5.  The  behavior  Is  similar  to  the  neg¬ 
ative-mass  behavior  that  one  would  expect  to  see  In  a 
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conventional  betatron  (Be  *  0).5  The  radial 

deflection  of  the  beam  seen  In  Fig.  5(a)  is 
accompanied  by  the  toroidal  bunching  In  Fig.  5(b). 
This  behavior  leads  us  to  compare  the  linear  growth 
rate  with  that  obtained  from  theory4  for  a  Bg  >  0 
beam.  We  find  that  the  theoretical  growth  rate  Is 
very  close:  1.63*10®  sec*1.  In  Fig.  6,  we  compare 
results  for  the  SL8  at  other  energies  with  the 
theory.  For  comparison,  results  for  the  same  beam  In 
a  modified  betatron  with  a  1  kG  toroidal  field  are 
shown.  In  Fig.  7,  the  same  curves  are  plotted  for  a 
1  kA  beam  with  a  *  1  cm.  For  these  parameters,  we 
have  just  one  simulation  result,  at  5.5  MeV,  and  the 
growth  rate  Is  again  seen  to  be  close  to  the  %  *  0 
result. 
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kA  SLB  beam.  Particles  at  all  e-positlons 
are  plotted  In  (a). 


We  must  note  that  the  results  In  Fig,  6  for  the 
case  8g  *  0  are  somewhat  artificial  when  y  <  18. 
In  this  regime,  weak  focusing  alone  Is  Insufficient 
to  hold  the  beam  together.  The  analytic  model  from 
which  the  growth  rates  are  obtained  uses  a  rigid  disk 
model  of  the  beam,  and  so  Ignores  the  force  balance 
required  within  the  beam.  It  Is  tempting  to  conclude 
from  the  results  In  Fig.  6  that  In  the  SLB,  the 
strong  focusing  provides  the  necessary  forces  the 
hold  the  beam  equilibrium  together,  but  that  the  neg¬ 
ative-mass  Instabilities  on  the  beam  behave  as  If  the 
toroidal  field  were  zero.  Further  numerical  and  ana¬ 
lytic  work  will  be  performed  to  check  this  conclus¬ 
ion. 
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Figure  6.  Growth  rates  of  t  *  20  mode  versus 

energy.  Simulation  results  for  SLB,  and 
analytic  results  for  conventional  (B8  = 
0)  and  modified  (Be  -  1  kG)  betatronsare 
shown. 
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Flqure  7.  Growth  rates  of  t  «  20  Instability  for  a  1 
kA  beam  with  1  cm  radius.  Curves  have 
same  meaning  as  In  Fig.  6. 
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APPENDIX  C 


STABILITY  OF  THE  SOLENOIDAL  LENS  BETATRON 


ABSTRACT 


The  solenoidal  lens  betatron  uses  solenoidal  focusing  to  increase  the 
current  which  can  be  injected  into  a  betatron.  A  dispersion  relation  for 
the  negative-mass  instability  in  this  device  is  derived  using  a  multiple- 
length-scale  method  to  average  over  the  nonuniform  toroidal  field.  The 
result  qualitatively  resembles  the  dispersion  relation  for  the  conventional 
betatron,  but  has  a  finite  transition  energy  and  suppressed  growth  rates. 
Both  effects  are  due  to  the  solenoidal  focusing.  It  is  shown  that  for  a 
space-charge-dominated  equilibrium,  finite  beam  radius  has  a  strong  stabi¬ 
lizing  influence  even  in  the  absence  of  any  energy  spread.  Results  of  3-D 
numerical  simulations  confirm  the  predicted  stability  at  low  toroidal  mode- 
numbers.  High  toroidal  mode-numbers  show  anomalously  large  growth  rates, 
but  saturate  in  a  relatively  benign  manner. 
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ILLUSTRATIONS 


Figure 


Conceptual  picture  of  the  sol enoi dal -lens  betatron.  In 
the  actual  device,  40  solenoids  are  placed  around  the 
racetrack 

The  coordinates  used  in  the  analytic  derivation  are 
cylindrical  (r,9,z),  and  toroidal  (p,<j>,0) 

Growth  rates  of  the  negative-mass  instability  obtained 
from  Eq.  (31)  for  the  parameters  in  Table  1 

An  equilibrium  in  which  there  is  little  transverse 
oscillation  of  the  particles  (as  in  a)  has  a  much 
greater  frequency  spread  than  an  equilibrium  in  which 
the  radius  is  determined  by  the  amplitude  of  betatron 
oscillations  (as  in  b) 

Branch  cut  in  the  complex  <u  -  tug  plane  to  make  the 
dispersion  relation,  Eq.  (31),  analytic.  The  arrows 
show  the  movement  of  the  roots  as  G  1  from  above 

Energy  history  of  the  i  *  15  mode  in  a  simulation  of  a 
100  A,  yq  =  3  beam,  showing  the  linear  growth  and  non¬ 
linear  saturation 

Particle  plots  (r-0,  z-9)  taken  during  the  nonlinear 
state  of  the  I  =  15  negative-mass  instability  on  a 
100  A,  yg  ■  3  beam.  The  dashed  lines  show  the  location 
of  the  toroidal ly  uniform  beam  at  the  start  of  the 
simulation 

Growth  rates  of  the  l  3  20  negative-mass  mode  on  a  10  kA 
beam  In  a  solenoldal-lens  betatron,  plotted  versus  beam 
energy.  Growth  rates  are  compared  to  those  for  a  con¬ 
ventional  betatron  (dashed  line)  and  those  obtained 
uslog  the  averaged-focusing  approximation  (solid  line) 

Particle  plots  (r-z,  r-0)  taken  during  the  nonlinear 
saturation  of  the  l  3  20  mode  on  a  10  kA  beam.  The  beam 
energy  is  yg  3  7  in  (a),  and  yg  ■  12  in  (b) 
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1.  INTRODUCTION 


Several  types  of  betatron-like  accelerators  for  high-current  charged 
particle  beams  are  being  studied  at  present  [1-5].  In  order  to  overcome 

•  the  space-charge  limit  on  the  current  at  low  energies,  these  devices  apply 
external  magnetic  fields  in  addition  to  those  of  the  conventional  betatron 

[6] .  At  the  Institute  for  Accelerator  and  Plasma  Beam  Technology  (IAPBT) 
of  the  University  of  New  Mexico,  a  device  is  under  construction  in  which 

•  the  additional  fields  take  the  form  of  solenoidal  lenses  [5],  as  shown  in 
Fig.  1.  Forty  solenoids  of  alternating  polarity  are  arranged  around  a 
racetrack-shaped  drift  tube.  Thus,  in  contrast  to  the  "modified"  betatron 
[1],  which  has  a  uniform  toroidal  magnetic  field,  the  toroidal  field  in  the 

•  IAPBT  betatron  alternates  in  direction  from  one  solenoid  to  the  next,  with 
a  magnetic  cusp  between  each  pair  of  solenoids.  The  advantages  of  this 
configuration  have  been  discussed  by  Humphries  et  al.  [5], 

!  In  this  paper,  we  consider  the  negative-mass  instability  [7]  in  the 

presence  of  solenoidal  lenses.  This  instability  has  been  shown  to  be 
important  for  high-current  electron  beams  in  the  modified  betatron  [8,  9], 
and  stellatron  [10],  potentially  causing  serious  disruption  during  the  time 

'  It  takes  to  accelerate  the  beam.  Derivation  of  a  negative-mass  instability 

dispersion  relation  for  the  solenoidal  lens  betatron  is  complicated  due  to 
the  lack  of  toroidal  symmetry.  To  deal  with  this  problem,  we  use  a 
multiple-length-scale  method  to  average  over  the  varying  solenoidal  fields. 

•  This  leads  to  a  dispersion  relation  like  that  for  the  conventional  betatron 

[7] ,  with  a  modified  transverse  focusing  term.  The  nature  of  the  beam 
equilibrium  in  the  solenoidal  lens  betatron  results  in  a  considerable 
spread  in  particle  circulation  frequencies  at  low  beam  energies.  We  evalu- 

•  ate  the  stabilizing  effect  this  has  on  the  negative  mass  instability.  The 
analytic  results  are  compared  to  three-dimensional  particle-in-cell  code 
simulations.  These  simulations  do  not  make  use  of  the  multiple-length- 
scale  approximation  and  so  can  be  run  at  short  wavelengths  where  this 

I  approximation  breaks  down. 


Conceptual  picture  of  the  solenoidal -lens  betatron.  In  the 
actual  device,  40  solenoids  are  placed  around  the  racetrack. 


This  paper  is  organized  as  follows.  In  Sec.  II,  we  obtain  the 
equations  for  the  equilibrium  and  small  amplitude  motion  of  the  beam.  In 
Sec.  Ill,  we  derive  a  dispersion  relation  for  the  negative-mass  instability 
including  the  effect  the  solenoidal  lenses.  In  Sec.  IV,  we  calculate  the 
effect  of  circulation  frequency  spread.  In  Sec.  V,  we  present  some  results 
of  particle  simulations  of  the  negative-mass  instability  in  a  solenoidal 
lens  betatron.  Our  conclusions  are  given  in  Sec.  VI. 

We  shall  use  a  system  of  normalized  units  which  are  convenient  for 

electron  beam  physics  problems.  A  plasma  frequency  Uq  is  defined  such  that 

c/ion  =  1  cm,  where  c  is  the  velocity  of  light.  Then  length  is  normalized 
u  2  2 

to  c/ojg ,  time  to  l/u>g,  velocities  to  c,  densities  to  u^m/4ue  ,  electric 

and  magnetic  fields  to  mcug/e,  where  e,  m  are  the  electronic  charge  and 

mass  respect i vely. 


2.  BEAM  EQUILIBRIUM 


For  the  purposes  of  this  paper,  we  will  assume  that  the  accelerator 
drift-tube  is  a  torus  rather  than  a  racetrack.  (The  IAPBT  device  has  a 
modular  design,  and  can  be  configured  either  as  a  racetrack  or  as  a  torus 
[5].)  The  coordinates  we  use  are  shown  in  Fig.  2.  The  reference  orbit  is 
that  of  a  particle  at  r  =  R,  z  =  0.  For  this  particle, 

R  =  V/Bz0  (i) 

where  B  _  is  the  value  of  the  vertical  field  at  r  =  R,  and  yn  is  the 
zu  2  wo  u 

matched  particle  energy  and  V  =  (1  -  1/yg)  '  .  If  we  displace  the  particle 

from  this  orbit,  then  its  equations  of  motion  are 


rVg  dp 

+  df  =  Er  +  Vz  -  VzB0 


VeVr  +  dP{ 
Y  —  +  dF 


=  E„  +  V  B  -  VB 
0  z  r  r  z 


dT  =  Ez  +  VrB0  -  V0Br 


where  the  components  of  the  vectors  V,  g,  E,  B,  represent  velocity, 
momentum,  electric  and  magnetic  fields,  respectively.  We  now  proceed  to 
linearize  these  equations  about  the  reference  orbit  in  the  paraxial 
approximation  [11].  For  a  particle  at  position  r  =  R  +  <Sr,  z  =  6z,  the 
applied  fields  are  [9] 

Bz  =  Bz0(l  -  n5r/R)  +  Bz$(e) 

8r  =  -Bz0  n«z/R  +  Brs(0) 

B0  "  B0s(9>  (3) 

where  n  is  the  external  field  index  n  =  -(R/Bzg)dBz/dr  and  subscript  s 
denotes  terms  due  to  the  solenoidal  magnets.  For  the  moment,  we  assume  the 
beam  is  in  a  region  of  uniform  toroidal  field  where  Brs>  Bzs  can  be 
neglected.  If  we  assume  a  Kapchinsky-Vladimi rsky  (KV)  [11]  beam  equi¬ 
librium  which  has  an  elliptical  cross-section  (we  will  show  that  this  is 
consistent  later)  then  the  self-fields  at  the  particle  can  be  calculated 
from  the  static  Maxwell's  equations: 

n  b 

E  =  ■  H  .  -  6r 

r  a  +  b 


a  +  b 


ngVb6r 


_  "p»*«r 
r  a  +  b 


where  a,  b  are  the  radii  in  the  r  and  z  directions,  respectively,  ng  is 
the  beam  density  and  superscript  s  denotes  "self-field".  Combining  these 
equations  with  the  linearizations  of  Eqs.  (2),  we  obtain 


T  T 


Sr  +  1  1  -  n  - 


+  6z 


0S 


6z  +  l  n 


6r 


fes 


*  o 


(5) 


2  3  2 

where  the  quantity  n$  =  2vR  /YgBgab  is  referred  to  as  the  self-field  index 
(v  is  Budker's  parameter).  For  a  self-consistent  KV  equilibrium,  the  radii 
a,  b  are  chosen  such  that 


2n$b  2nsa 

a  +  b  p  a  +  b 


(6) 


To  treat  the  effect  of  the  magnetic  cusps,  it  is  convenient  to  change  the 

coordinate  system  to  the  toroidal  coordinates  shown  in  Fig.  2  and  introduce 

a  vector  potential  A^Je)  which  describes  the  solenoidal  field.  Equations 

$s 

(5)  are  replaced  by 


•2  2  •  0s  „ 

p  -  p<|>  +  w  p  +  p$  - —  =  0  , 


tp<4  *  Vl  =  0  • 


2  (l  n  2"sb  )  2 

%  =  I1  -  n  ”  T"+~b  /  0  • 


Equation  (7b)  can  be  integrated  to  obtain 


where 


,  2*  pAfs  2  (•  l  Bes' 

L  =  P  $  +  - -  P  1  $  -  -rr  - 


’0 


const. 


(7a) 

(7b) 


(8) 


Since  L  is  a  constant,  its  value  is  determined  by  the  Initial  conditions. 
Experimentally,  the  beam  will  be  generated  using  a  cathode  shielded  from 
magnetic  fields  [5],  so  that  1=0.  Substituting  this  into  Eq.  (7a),  we 
obtain 
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P  *  0  . 


(9) 
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In  the  limit  that  the  cusps  become  infinitely  sharp,  B0S  is  a  constant  over 
the  particle  orbit,  and  Eq.  (9)  has  a  stationary  solution  for 


B2 

1  Bes  .  2 

T  ~2~  “p 

Yo 


=  0 


(10) 


This  is  the  condition  for  a  space-charge  dominated  equilibrium,  l.e.,  one 
where  the  emittance  is  negligible.  To  be  consistent  with  the  assumption  of 
a  fixed  equilibrium  profile,  we  require  that  the  "phase-advance  per  sole¬ 
noid”,  i.e.,  the  poloidal  angle  through  which  the  particles  rotate  on  gong 
through  one  solenoid,  be  small: 


u0 


l  Bes  s 


« 


l 


(ID 


where  S  is  the  length  of  a  solenoid.  The  case  n  *  1/2  is  a  special  one, 
since  a  ■  b,  and  poloidal  rotation  does  not  affect  the  beam  profile.  In 
this  case,  Eq.  (11)  can  be  relaxed.  Struckmeler  and  Reiser  [12]  have  shown 
that  uq  <  w/2  is  required  for  envelope  stability. 


Physically,  the  equilibrium  we  have  just  constructed  behaves  as 
follows.  In  a  given  solenoid,  the  beam  particles  rotate  pololdally  at  the 
Larmor  frequency,  B0s/2yo<  On  passing  through  the  magnetic  cusp  Into  the 
next  solenoid,  the  poloidal  rotation  changes  direction.  The  shape  of  the 
beam  Is  elliptical.  In  general,  with  radii  determined  from  Eq.  (6),  and  has 
a  uniform  density  determined  from  Eq.  (10).  If  the  magnetic  cusps  are  not 
Infinitely  sharp,  then  an  equilibrium  can  stIH  be  constructed  [12],  but 
the  radii  a,  b  vary  periodically  with  toroidal  angle. 
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3.  NEGATIVE-MASS  INSTABILITY 


In  order  to  analyse  the  negative-mass  instability,  we  adopt  a 
simplified  model  of  the  beam  dynamics.  We  assume  that  the  beam  acts  like  a 
rigid  body  in  the  transverse  direction,  so  that  only  the  transverse  motion 
of  the  beam  centroid  has  to  be  calculated,  rather  than  the  transverse 
motion  of  each  beam  particle.  Comparisons  between  a  rigid-beam  model  and 
three-dimensional  particle  simulations  for  the  case  of  the  modified  beta¬ 
tron  [8]  have  produced  good  agreement.  The  physical  basis  for  the  model 
rests  on  the  fact  that  when  the  beam  minor  radius  is  much  smaller  than  the 
drift  tube  minor  radius,  the  perturbed  toroidal  forces  are  the  same  for  all 
particles  in  a  given  transverse  slice  of  the  beam,  and  the  perturbed  trans¬ 
verse  forces  are  linearly  proportional  to  the  transverse  displacement  of 
each  particle  in  the  slice.  In  the  toroidal  direction,  the  model  places  no 
restriction  on  the  beam  motion.  The  transverse  motion  of  the  beam  gener¬ 
ates  perturbed  dipole  fields,  while  toroidal  bunching  generates  perturbed 
monopole  fields. 

To  obtain  the  equations  for  the  beam  centroid,  we  start  from  Eqs.  (2). 
The  field  quantities  which  we  Insert  into  these  equations  are  evaluated  at 
the  center  of  the  beam.  It  has  been  shown  [8]  that  for  high-current  beams, 
toroidal  corrections  to  the  continuity  equation  and  the  field  equations 
(l.e.,  terms  of  order  d/R  and  higher  where  d  Is  the  wall  radius)  must  be 
retained  in  the  calculation,  particularly  for  high  toroidal  mode  numbers. 

A  model  which  uses  an  exact,  though  cumbersome,  solution  to  these  equations 
has  been  derived  [13],  The  exact  model  shows  that  the  cylindrical  approxi¬ 
mation  for  the  field  equations  gives  reasonably  good  results  for  low  toroi¬ 
dal  mode-numbers  If  a  certain  replacement  is  made  [Eq.  (30)].  Here,  we 
will  make  use  of  this  simplification  in  order  to  avoid  the  complexities  of 
toroidal  corrections. 

For  low  toroidal  mode-numbers,  the  long-wavelength,  low-frequency 
2  2  2  2 

approximation  u  ,  t  /R  «  ?f  (u,  I  denote  the  mode  frequency  and 
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wave-number,  respectively,  and  denote  the  transverse  Laplacian)  can  be 
used  in  solving  the  field  equations.  Neglecting  toroidal  corrections  and 
assuming  a  circular  beam  cross-section,  the  solutions  for  the  perturbed 
transverse  self-fields  at  the  perturbed  beam  centroid  position  are 


c  n0  a2 
AEr  =  Ar 


"0  a2 

lEz  -T  Ja 


s  n0V  a2 
AB^  =  Ar 

r  2  d2 


c  n0V  2 

4BI  *-T7“ 


where  A  denotes  perturbed  quantities,  and  Ar,  az  denote  the  transverse  dis¬ 
placement  amplitudes  of  the  beam  centroid.  These  fields  come  from  the 
dipole  components  of  aEj_  and  AB^.  The  monopole  components  of  the  trans¬ 
verse  fields  do  not  couple  to  the  centroid  motion.  The  monopole  component 
of  aEq  does  enter,  however,  and  can  be  obtained  using  the  integral  form  of 
Faraday's  Law, 


f  1  • « -  -  It  /  » •  « 

lEe  '  ‘lie  /  4E>  -  It  /  4B$d» 


where  the  superscript  zero  denotes  the  monopole  component.  The  sources  for 
AE®,  AB^  are  the  monopole  components  of  An  and  AJ0,  respectively.  The 
latter  are  related  to  the  perturbed  beam  displacement  through  the  con¬ 
tinuity  equation  [8]: 
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4"  '  -  I?  Vz  -  F  If  r"o4r  -  Te  V9 


From  this  and  aJq  =  nQRA9  +  AnV  where  the  dot  denotes  the  total  time 


derivative,  we  obtain  the  expressions 


4"°  -  - 1?  v9 


Al0  3 

Ajz  =  at  V6 


(15) 


Consistent  with  our  use  of  the  cylindrical  approximation,  we  have  neglected 
toroidal  terms  in  Eq.  (15).  Substitution  into  Eq.  (14)  yields 


4EMv2  (lt2Mf)(^ 


(16) 


The  applied  fields  are  given  by  Eq.  [3],  with  5r,  Sz  replaced  by  Ar,  az. 
We  can  now  write  down  the  linearized  equations  of  motion  for  the  beam 
centroid: 


Ar 


r  +  4Ar  +  4z  ^  =  y2 u  aV< 

Y0  T0 


'0  0°  0 


2  •  B8s  AVeBrs 

Az  +  wJAZ  -  Ar  + 


Y0  Y0 


=  0 


4«.-t 

- 


(17) 


where 


and  we  have  taken  n  *  1/2  to  be  consistent  with  a  =  b  [Eq.  (6)].  To  solve 
Eq.  (17),  we  perform  a  multiple-length-scale  analysis  to  average  over  the 
solenoidal  lenses.  This  Involves  an  ordering  scheme  where  the  wavelength 
of  the  beam  mode  \  is  assumed  to  be  much  longer  than  the  period  of  the 
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alternating  solenoidal  fields,  i.e.,  S/x  =  0(e),  where  e  is  a  smallness 
parameter.  Consistent  with  this,  we  assume  3/3t,  «jg,  wj_  are  of  order  e. 
In  addition,  we  assume,  as  in  the  equilibrium  equations,  that  Ug  <  1 
[Eq.  (11)].  It  is  convenient  to  perform  the  averaging  in  toroidal  coordi¬ 
nates  (Fig.  2),  in  which  Eqs.  (17)  have  the  form 


•2.2,;  ®0s  2  ... 

P  -  P$  +  <*>JP  +  y -  “  Y0“0AV9  C0Sl*’ 


(18a) 


p  dt  p  [* 


1  B9s 


=  ^0“oiV9  Sin  ♦ 


(18b) 


C  +  io,,Ar  -  C 


2  2 

~~1 - ? 

3S  3t 


5  =  0 


(18c) 


where  c  =  RA0,  s  =  Re,  and  C  =  l/4(nga^)(l  +  2*n  d/a)/Yg.  All  quantities 
are  expanded  in  powers  of  e  in  the  following  manner. 


p(s,t)  *  pg(sQ,s^,t)  +  epj(sQ,Sj,t)  + 


(19a) 


♦(s.t)  ■  ♦0(s0,s1,t)  +  •♦1(s0.s1.t)  + 


(19b) 


C(s,t)  -  CgCsg.Sj.t)  +  eCjCsg.Sj.t)  + 


(19c) 


where  sQ,  Sj  denote  the  short  and  long  length-scales,  respectively.  The 
essence  of  the  multiple-scale  method  is  that  sQ  and  s^  are  treated  as  inde¬ 
pendent  variables.  The  total  time  derivative  d/dt  Is  expanded  as 
V(3/3s0)  +  e(d/dtj)  where  d/dt^  *  3/ at  +  Va/3s^.  To  zero  order  in  e, 

Eqs.  (18)  become 
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From  these,  we  conclude  that  pQ,  $q,  Cq  do  not  depend  on  Sq.  Using  this 
information,  we  obtain  to  order  e. 
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These  equations  imply  that  p^,  ^  also  do  not  depend  on  Sq.  For  <^,  we 
obtain 


/p(t,s^)  3$q\  1  (  ®es 

*  U  aT" )  S0  +  ?  J  TjJ"  ds0  +  G^*sl  ) 


(22) 


where  Ftt.Sj),  G(t.Sj)  are  arbitrary  constants.  The  secular  dependence  on 
Sq  can  be  removed  by  choosing 


2  3*0 

F  *  po  ir 


(23) 


This  choice  implies  that  d^j/asQ  depends  only  on  Sq.  Proceeding  to  order 


e  ,  we  obtain 
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We  now  average  these  equations  over  one  period  of  the  magnetic  lenses. 
Subtracting  the  averaged  equations  for  Eqs.  (24b)  and  (24c)  from  the 
unaveraged  ones,  we  obtain 


(25a) 


(25b) 


where  the  brackets  <>  denote  the  average.  Solving  these  equations,  and 

2  2  2  2 

removing  secularltles  In  the  usual  way,  we  find  a’^/^Sg  *  3  ^/asg  *  0. 
From  Eq.  (24b),  we  note  that  If  we  Impose  the  initial  conditions  *  0, 

340/3t  ■  0  at  t  *  0,  then  this  equation  has  the  unique  solution  4g(t)  *  0. 
By  Imposing  these  Initial  conditions,  and  thereby  restricting  the  class  of 
Initial  conditions  for  which  the  analysis  is  valid,  we  can  set  4>0  -  0  In 
Eqs.  (24).  From  Eqs.  (22)  and  (24a),  we  then  obtain 
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By  subtracting  this  equation  from  its  average,  it  is  again  straightforward 
to  show  that  <32p2/3Sq>  =  0  (but  32p2/3s2  *  0).  We  thereby  obtain  the 
equations  for  the  averaged  beam  centroid  motion: 


V.  As 
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(27a) 
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(27b) 


(27c) 


Note  that  Eqs.  (27a)  and  (27b)  have  the  form  of  the  equations  one  would 

obtain  if  no  solenoidal  magnetic  lenses  were  present.  Thus,  to  the  order 

we  have  solved  the  original  equations,  the  net  effect  of  the  solenoidal 

lenses  is  to  increase  the  transverse  focusing  on  the  beam  by  an  amount 
2  2 

proportional  to  <B0s/Yq>,  i.e.. 


4*“i  +  ni 


To  obtain  -a  dispersion  relation  from  Eqs,  (27),  we  assume  that  the 

perturbed  quantities  vary  as  exp(1Z0  -  lut).  In  Eq.  (27b),  we  obtain  the 
2  2  2 

factor  ft  /R  -  «  .  However,  as  shown  in  Ref.  [13],  when  toroidal  correc¬ 
tions  to  the  field  equations  and  the  continuity  equation  are  kept,  we 

2  2  2 

Instead  obtain  the  factor  c^ft  /R  -  ,  where  and  a,  are  frequency- 

and  mode-number-dependent  expressions  approximately  equal  to  unity.  In 
general,  we  can  approximate  u  -  ft^g  in  evaluating  this  factor. 
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Small  differences  between  and  can  strongly  affect  the  magnitude  and 
even  the  sign  of  -  c^V^.  These  effects  become  increasingly  evident  as 
the  mode  number  1  increases.  However,  the  averaging  procedure  we  have  used 
is  only  valid  for  low  mode  numbers,  and  for  our  purposes,  it  will  be  suffi¬ 
cient  to  set  o1  =  a2  =  1  in  Eq.  (29) 
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This  approximation  was  made  on  heuristic  grounds  by  Landau  and  Neil  [7], 

The  main  error  introduced  is  in  the  asymptotic  fall-off  of  growth  rate  with 

Yn  [13].  The  important  point  here  is  that  much  larger  errors  are  intro- 

u  2  2  2 
duced  by  keeping  the  unmodified  factor  l  /R  -  u»  in  the  dispersion  rela¬ 
tion.  Fictitious  cutoffs  in  the  growth  rate  as  a  function  of  energy  and 
current  are  predicted  [14]. 


Using  Eq.  (28),  the  dispersion  relation  obtained  is: 


1  = 


1 


(31) 


where  n  a  u>  -  tu>g. 
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4.  SOLUTIONS  TO  THE  DISPERSION  RELATION 


2 

By  inspection  of  Eq.  (31),  we  see  that  there  are  resonances  at  ft  =  0 
and  ft2  «  ft£,  corresponding  to  the  longitudinal  and  transverse  modes  of 
oscillation  of  the  beam.  These  modes  are  coupled  through  the  negative-mass 
effect.  For  the  low  i-numbers  we  are  considering,  the  coupling  is  non¬ 
resonant,  i.e.,  the  resonances  do  not  overlap.  Thus,  for  the  mode  near 
ft  »  0,  we  can  set  ft2  «  ft|_,  to  get 


_2  .  t2C  /l  “0 


When  the  right-hand  side  is  negative,  we  obtain  the  negative-mass  insta¬ 
bility.  For  B0s  =  0,  Eq.  (32)  becomes  the  expression  of  Landau  and  Neil 
[7]  for  a  monoenergetic  beam  in  a  conventional  betatron.  The  condition  for 
instability  is 


Thus,  the  strong-focusing  effect  of  the  solenoids  introduces  a  finite 
negative-mass  transition  energy  Ytr  below  which  the  beam  is  stable.  This 
contrasts  with  the  conventional  betatron  which  is  unstable  at  all  energies. 
Further,  we  note  from  Eq.  (32)  that  when  <b|  /yq>  >  1/2,  the  growth  rate 
scales  as  1/ j BQs j ,  so  that  the  solenoidal  focusing  strongly  suppresses 
growth  of  the  testability. 


As  the  beam  is  accelerated,  Yg  Increases.  We  see  from  Eq.  (28)  that 
the  effect  of  the  solenoidal  focusing  rapidly  decreases  unless  B0S  is 
increased  also.  By  ramping  B0g  In  synchronism  with  Yq.  the  suppression  of 
growth  rates  relative  to  the  conventional  betatron  is  maintained.  An  addi¬ 
tional  benefit  Is  that  the  crossing  of  single-particle  resonances  is 
avoided.  The  obvious  drawback  Is  that  additional  energy  is  required  to 
drive  the  solenoids. 


In  Fig.  3,  we  show  growth  rates  obtained  from  Eq.  (29)  for  parameters 
typical  of  the  IAPBT  betatron  (Table  1).  The  curves  show  the  contrast 
between  a  case  where  the  solenoidal  field  is  kept  fixed  at  its  injection 
value,  and  one  where  it  is  ramped  linearly  with  Yq.  Note  that  the  growth 
rates  in  both  cases  are  substantial  for  a  device  an  acceleration  time  of 
several  tens  of  microseconds. 


TABLE  1.  SOLENOIDAL  LENS  BETATRON  PARAMETERS 


» 

Torus  Major  Radius  1  m 

Torus  Minor  Radius  4.5  cm 

I  Beam  Major  Radius  1  m 

Beam  Minor  Radius  (at  injection)  -2  cm 

Injected  Beam  Energy  300  kV 

(  Injected  Beam  Current  “100  A 

Solenoidal  Field  Amplitude  (at  injection)  -220  Gauss 

Number  of  Solenoids  30 


There  is  a  qualitative  difference  between  Eq.  (31)  and  the  the  corre¬ 
sponding  dispersion  relation  for  the  modified  betatron  [1,  9],  which  has  a 
uniform  toroidal  field.  In  Eq.  (31)  only  the  radial  and  toroidal  modes  of 
the  beam  are  coupled  (cf.  Fig.  7).  In  the  modified  betatron,  the  radial 
beam  motion  is -coupled  to  the  vertical  motion  (z-direction  in  Fig.  2)  as 
well  as  to  the  toroidal  motion. 
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Figure  3.  Growth  rates  of  the  negative-mass  instability  obtained  from 
Eq.  (31)  for  the  parameters  in  Table  1. 
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5.  STABILIZATION  DUE  TO  FINITE  BEAM  RADIUS 


In  the  derivation  of  Eq.  (31),  we  assumed  that  all  particles  rotate 
about  the  major  axis  of  the  torus  at  the  same  frequency,  namely,  =  V/R. 
For  a  finite  radius  beam,  however,  there  is  a  spread  in  path  lengths  around 
the  torus.  Thus,  even  if  V  =  c  for  all  particles,  there  is  a  spread  in 
rotation  frequencies.  Taking  a  uniform  density  equilibrium  of  the  type 
discussed  in  Sec.  II,  with  circular  cross-section  and  small  phase  advance 
•  per  solenoid  (Eq.  11),  it  is  strai ghtforward  to  show  that  the  distribution 

function  for  the  particle  rotation  frequencies  is 

f(9)  =  A  [a2*2/R2  -  (9  -  .c)2]1/2  (34) 

I 

where  A  is  a  normal i zation  constant  and  9  is  the  particle  rotation  fre¬ 
quency.  The  width  of  the  distribution  function  is  awg/R.  This  contrasts 

with  the  conventional  and  modified  betatrons,  where  the  betatron  oscilla- 

2 

Lions  of  the  particles  result  in  a  width  which  is  only  of  order  (a/R)  for 
a  monoenergetic  beam  (see  Fig.  4).  For  these  devices,  it  is  generally 


I  Figure  4.  An  equilibrium  in  which  there  is  little  transverse  oscillation 

of  the  particles  (as  in  a)  has  a  much  greater  frequency  spread 
than  an  equilibrium  in  which  the  radius  is  determined  by  the 
amplitude  of  betatron  oscillations  (as  in  b). 
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necessary  to  have  an  energy  spread  on  the  beam  to  obtain  a  significant 
frequency  spread  [7,  9]. 


i 


I 


» 


» 


I 


In  order  to  derive  a  dispersion  relation  taking  frequency  spread  into 
account,  a  kinetic  treatment  is  necessary.  Rather  than  redoing  the 
analysis  of  Sec.  Ill  in  the  framework  of  the  Vlasov  equation,  we  will  make 
a  plausibility  argument  for  adapting  the  kinetic  dispersion  relation  for 
the  conventional  betatron  derived  by  Landau  and  Neil  [7].  As  observed  in 
Sec.  IV,  Eq.  (31)  can  be  obtained  from  the  cold  beam  result  of  Ref.  [7] 

when  the  replacement  in  Eq.  (28)  is  made.  The  longitudinal  motion  of  the 

2 

particles  is  affected  by  the  solenoidal  field  only  to  order  (uQa/S)  «  l. 

Therefore,  we  expect  the  influence  of  the  solenoids  on  the  resonance  at 

2 

fl=0,  which  is  associated  with  the  instability,  to  be  minimal.  In  addi¬ 
tion,  as  noted  at  the  beginning  of  Sec.  4,  there  is  no  resonant  coupling  to 
the  transverse  beam  modes  for  low  ^-numbers.  On  this  basis,  we  modify  the 
kinetic  equation  of  Ref.  [7]  according  to  Eq.  (28),  we  obtain  the  disper¬ 
sion  relation: 
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where 


(35) 


2 

p  =  -yujqR  +  p  is  the  canonical  angular  momentum,  and  f  is  normalized  such 
that  / fdp  =  ng.  The  relation  between  p  and  9  is  9  =  ui^  -  kp.  Converting 
Eq.  (35)  into  an  integral  over  9,  and  inserting  the  distribution  function 
in  Eq.  (34),  we  obtain  the  dispersion  relation 
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(36) 


I 


C-24 


where 


G 


Equation  (36)  can  be  made  analytic  In  the  complex  w  plane  by  Introducing 
the  branch  cut  shown  in  Fig.  5.  We  find  that  for  1  <  G  <  the  beam  is 
unstable  with  growth  rate 


r 


lu0a  G  -  1 

"T~”  mr: r 


(37) 


The  zero-frequency-spread  growth  rate  Tg  In  Eq.  (32)  can  be  recovered  by 
letting  a  ♦  0.  As  G  approaches  1  from  above,  the  roots  approach  the  branch 
cut,  as  shown  In  Fig.  5.  For  G  <  1,  the  negative-mass  mode  disappears,  and 
the  beam  Is  stable.  The  stability  criterion  can  be  rewritten 


l“0d 


(38) 


This  criterion  has  a  physical  Interpretation  In  terms  of  the  rotation 
frequency  spread  of  the  particles.  Consider  two  particles,  one  at  r  ■  R, 
and  one  at  r  *  R  t  a,  which  are  Initially  at  the  same  toro’dal  angle.  Then 
at  marginal  stability,  these  particles  will  become  separated  by  1 / 1  radians 
in  a  time  l/Tg.  This  smearing  effect  prevents  clumping  of  the  beam. 


We  emphasize  that  the  stability  criterion  just  derived  is  only  valid 
for  a  cold,  space-charge  dominated  equilibrium  of  the  type  described  in 
Sec.  II.  To  apply  Eq.  (38)  to  a  hot  beam,  one  must  calculate  an  effective 
radius  which  excludes  the  contribution  of  transverse  betatron  oscillations 
to  the  radius.  In  the  extreme  case  where  the  finite  beam  radius  is  com¬ 
pletely  due  to  betatron  oecl 1 lat ions ,  as  in  Fig.  4b,  the  frequency  spread 
1$  negligible  f or  typical  parameters. 
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/-.plying  Eq.  (38)  to  the  cases  in  Fig.  3,  we  find  that  a  beam  radius 
of  1.2  cm  is  sufficient  to  stabilize  the  peak  growth  rate.  As  the  beam  is 
accelerated,  its  radius  shrinks,  thereby  decreasing  the  width  of  the 
frequency  distribution.  From  Eq.  (10)  we  find 


-i 


(39) 


so  that  the  radius  shrinks  as  yq  if  <Bftc/yn>  is  held  fixed.  If  B„_  is 


.-3/2  „  ,q2  y  2 

es/Yo  "  -  •  . . .  “es 

held  fixed,  then  the  scaling  is  more  complicated  at  low  energy,  but  quickly 
3/2 

goes  over  to  a  Yn  scaling  as  Yn  increases.  From  Eq.  (32),  we  see  that 
u  u  3/2 

the  growth  rate  rQ  also  shrinks  as  yq  .  Therefore,  if  the  beam  is  stable 
just  above  the  transition  energy,  then  it  will  tend  to  remain  so  as  the 
acceleration  proceeds.  In  addition,  both  sides  of  Eq.  (38)  scale  linearly 
with  t,  so  that  if  one  mode  is  stable  then  all  modes  are.  As  we  shall  see 
in  Sec.  6.1,  this  scaling  breaks  down  for  large  £-numbers. 


C-27 


6.  COMPARISON  WITH  NUMERICAL  SIMULATIONS 


In  order  to  check  the  analytic  theory,  and  to  study  the  nonlinear 
effects  of  the  instability,  we  carried  out  a  limited  number  of  simulations 
'  using  the  three-dimensional  particle-in-cell  code  IVORY,  IVORY  has  been 

used  successfully  to  model  the  negative-mass  instability  in  the  modified 
betatron  [8,  13],  The  code  is  fully  electromagnetic  and  advances  the  com¬ 
plete  Lorentz  force  equations  for  the  particles.  A  two-dimensional  grid  is 
»  used  to  represent  the  transverse  plane  of  the  beam,  while  fields  in  the 

toroidal  direction  are  represented  by  a  Fourier  sum.  This  allows  us  to 
compare  predictions  for  individual  4-numbers  directly  with  theory.  For 
economy,  we  generally  keep  < 3  toroidal  modes  in  the  simulations.  To  model 

I  the  solenoidal  lens  betatron,  a  sinusoidally  varying  toroidal  field  is 

2  2  2  2 

used.  We  use  Eq.  (10),  with  BQs/yq  is  replaced  by  <Bgs/YQ>.  to  obtain 
suitable  equilibria  to  initialize  the  simulations.  We  find  that  small- 
amplitude  envelope  oscillations  are  excited  due  to  the  fact  that  the  cusps 
*  are  not  infinitely  sharp  [15], 

6.1  IAPBT  PARAMETERS 

Three  simulations  were  performed  for  the  for  the  parameters  In 
Table  1.  We  assumed  30  solenoids  around  the  torus  since  10  of  the  40  sole¬ 
noids  In  the  actual  device  are  on  straight  sections  which  we  do  not  simu¬ 
late  here.  First,  we  looked  at  the  t  *  5  mode  at  yQ  ■  3.  The  beam  radius 
»  was  chosen  to  be  1.75  cm,  for  which  the  matched  amplitude  of  the  solenoidal 

field  Is  134  gauss.  The  simulation  was  run  for  230  ns,  during  which  no 
measurable  growth  was  observed.  The  theoretically  predicted  growth  rate  is 
22.8  us‘l,  which  would  have  produced  about  5.3  e-foldings  In  230  ns.  (This 
I  growth  rate  Is  higher  than  those  shown  In  Fig.  3  because  of  the  lower 

solenoidal  field.)  We  attribute  the  beam  stability  to  the  spread  In 
circulation  frequencies.  The  right-hand  side  of  Eq.  (38)  Is  26.2  us*1, 
which  Is  marginally  sufficient  for  stability. 
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The  averaging  procedure  of  Sec.  4  breaks  down  as  the  mode-number  t 
approaches  the  mode-number  associated  with  the  solenoldal  field.  For 
our  case,  *  IS.  Since  no  averaging  Is  used  In  the  simulation  code, 
however,  such  modes  can  be  looked  at  numerically.  For  the  same  parameters 

as  the  i  *  5  simulation  above,  we  performed  a  simulation  of  the  t  ■  1$  «  15 
mode.  We  observed  a  strong  Instability,  with  a  growth  rate  of  66  us-*,  as 
shown  In  Fig.  6.  A  particle  plot  from  the  nonlinear  stage  of  the  insta¬ 
bility  (Fig.  7)  shows  typical  negatl ve-mass  Instability  behavior,  l.e., 
radial  kinking  accompanied  by  toroidal  clumping.  The  Instability  has 
almost  saturated  at  the  point  reached  In  Fig.  7,  but  1/6  of  the  current  has 
been  scraped  off  onto  the  wall.  If  the  wall  radius  ha  1  been  larger,  the 
Instability  might  have  saturated  without  current  loss. 

If  we  naively  apply  the  analytic  dispersion  relation  to  the  t  *  15 
mode,  we  obtain  a  growth  rate  of  66  us"1,  the  same  as  the  numerical  result. 
However,  Eq.  (38)  predicts  stability  due  to  frequency  spread.  Thus,  It 
appears  that  the  zero-frequency-spread  growth  rate  for  1  ■  15  Is  signifi¬ 
cantly  larger  than  that  predicted  by  analytic  theory.  On  the  other  hand, 


Figure  6.  Energy  history  of  the  t  *  15  mode  In  a  simulation  of  a  100  A, 
TO  “  3  beam,  showing  the  linear  growth  and  nonlinear 
saturation. 
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the  analytic  prediction  of  a  negative-mass  transition  energy  appears  to  be 
valid  even  for  i  *  15.  We  simulated  this  mode  for  Yg  =  1.6,  which  corre¬ 
sponds  to  the  injection  energy  of  the  IAPBT  betatron,  and  observe  no  insta¬ 
bility.  The  analytic  theory  predicts  stability  up  to  at  least  y0  *  2.2 
(see  Fig.  3). 

6.2  HIGH-CURRENT  BEAM  STABILITY 

If  experiments  at  100  A  are  successful,  then  higher-current  experi¬ 
ments  will  be  undertaken.  We  have  performed  three  simulations  of  10  kA 
beams  to  Investigate  beam  stability  at  high  current.  We  assumed  20  sole¬ 
noids  around  a  torus  with  aim  major  radius,  and  simulated  the  mode  l  *  2 
*  20  (this  mode  was  chosen  for  reasons  of  computational  economy).  In 
Fig.  8,  we  compare  the  simulation  results  with  (a)  growth  rates  obtained 
for  a  conventional  betatron  (no  toroidal  field)  and  (b)  growth  rates 
obtained  by  using  the  averaged  approximation  for  the  solenoidal  lenses. 

Because  of  the  high  t-number,  the  curves  In  Fig.  8  were  obtained  by  insert¬ 
ing  Eq.  (28)  Into  the  exact  analytic  dispersion  relation  of  Ref.  (13) 
rather  than  using  Eq.  (31).  In  doing  the  simulations,  we  set  the  beam 
radius  to  2  cm  In  order  to  have  good  resolution  of  t  e  beam  on  the  simula¬ 
tion  mesh.  This  means  that  the  solenoidal  field  necessary  to  get  a  matched 
beam  decreases  as  the  beam  energy  Increases  (see  Eq.  39).  Thus  the  simu¬ 
lation  growth  rates  shown  In  Fig.  8  do  not  represent  those  experienced  by 
an  accelerating  beam  where  the  solenoidal  field  Is  kept  constant,  or 
Increased. 

As  In  the  100  A  simulations,  we  find  that  the  hlgh-t  growth  rates  are 
anomalously  large  when  frequency  spread  effects  are  taken  Into  account. 

Nevertheless,  the  simulations  for  y0  ■  7,  and  y0  *  12  both  saturate  with  no 
loss  of  current.  Plots  from  the  nonlinear  stage  of  these  simulations  are 
shown  In  Fig.  9.  There  Is  considerable  churning  of  the  beam  after  these 
plots  were  taken,  and  this  stage  of  the  Instability  Is  probably  not  well 
modeled  with  a  few  Fourier  modes.  The  simulation  at  yq  *  17  was  not  run 
long  enough  to  see  saturation. 

y 

i 

I 


I 


C-31 


8.80  -*».*» 0  0.00  4.40  8.80 

i  (cm) 


91 


96  100  104  109 

r  (cm) 


R-9  22 


Figure  9.  Particle  plots  (r-t,  r-e)  taken  during  the  nonlinear  saturation 
of  the  t  *  20  mode  on  a  10  kA  beam.  The  beam  energy  Is  yg  •  7 
In  (a),  and  yn  •  12  In  (b). 


7.  SUMMARY  AND  DISCUSSION 


We  have  presented  an  analytic  model  for  the  negative-mass  instability 
in  the  solenoidal -lens  betatron.  Restricting  the  analysis  to  low  toroidal 
mode-numbers,  the  effect  of  the  solenoidal  lenses  are  averaged  over  using  a 
multiple-length-scale  method.  The  equations  obtained  are  those  of  the  con¬ 
ventional  negative-mass  dispersion  relation  with  a  modified  transverse 
focusing  force.  The  additional  transverse  focusing  introduces  a  finite 
negative-mass  transition  energy  below  which  the  beam  is  stable.  Above  this 
energy,  growth  rates  are  suppressed  relative  to  those  of  the  conventional 
betatron. 

The  averaged  transverse  focusing  term  is  used  to  obtain  a  kinetic 
dispersion  relation  for  the  solenoidal  lens  betatron.  We  find  that  for  a 
space-charge-dominated  equilibrium,  it  is  relatively  easy  to  stabilize  tne 
Instability  with  the  natural  frequency  spread  due  to  finite  beam  radius. 
This  contrasts  with  other  types  of  betatrons  where  an  energy  spread  is 
needed  to  produce  a  significant  frequency  spread. 

Results  of  three-dimensional  PIC  code  simulations  of  the  IAPBT  beta¬ 
tron  are  in  agreement  with  the  predicted  stability  for  low  t-numbers  where 
the  averaging  procedure  Is  applicable.  For  a  mode-number  equal  to  the 
mode-number  of  the  solenoidal  lenses,  on  the  other  hand,  we  find  that  the 
growth  rate  is  anomalously  large.  In  the  nonlinear  regime,  the  instability 
exhibits  a  classical  wave-breaking  saturation  mechanism  leading  to  some 
loss  of  current. 

The  prediction  of  a  finite  negative-mass  transition  energy  Is  borne 
out  by  the  simulation  results  even  for  large  I -numbers.  Thus,  the  IAPBT 
betatron  Is  predicted  to  be  stable  at  its  Injection  energy  even  In  the 
absence  of  particle  frequency  spread. 
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High-current  simulations  at  high  mode-number  show  large  growth  rates 
comparable  to  what  one  would  obtain  in  the  absence  of  any  toroidal  field 
(i.e.,  a  conventional  betatron).  However,  no  current  loss  is  observed 
during  the  nonlinear  saturation  of  the  Instability. 
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